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The range of fission fragments in different gases has been studied by means of the cloud- 
chamber method. The reduced range was found to increase from xenon to helium, and the 
range in helium was found to be about 18 percent larger than in argon, which is in agreement 
with previous determinations. The reduced range in hydrogen was found to be smaller than in 
other gases, pointing to a comparatively larger electronic stopping effect. The reduced ranges 
of fragments passing through hydrogen and deuterium, respectively, have been compared and 
the range in deuterium was found to be about 6 percent larger than that in hydrogen. Since the 
electronic stopping effect is identical in both gases, the difference in the ranges must be 
ascribed to the nuclear stopping effect. The spread of the range within the groups was discussed 
in relation to the energy measurements by Flammersfeld, Jensen, and Gentner, using Bohr’s 
theory of the stopping effect, and a reasonable agreement was obtained. The previous con- 
clusions concerning the velocity-range relations of fission fragments were confirmed by new, 
more direct evidence, using close collisions for a determination of the velocity or mass of the 


fragments. 


AS reported in a previous paper,! the ratio 
between the range of fission fragments and 
the range of a-rays in helium was found to be 
about 20 percent larger than the corresponding 
ratio in argon. As pointed out by Bohr,’ this 
difference is just in agreement with the different 
shapes of the velocity-range curves in light and 
in heavy gases to be expected from theoretical 
considerations concerning the electronic and 
nuclear stopping effects. A closer examination® 
of the range-velocity relation of fission fragments 
in helium and a comparison with the corre- 
sponding relation in argon has, within the experi- 
mental errors, confirmed the theoretical con- 


siderations and indicates that the rate of. 
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velocity loss is similar in the two gases in the first 
part of the range, where the stopping is caused 
mainly by electronic encounters, whereas the end 
part of the velocity-range curve in helium is 
comparatively longer than that in argon, cor- 
responding to the effect of the smaller energy 
loss by nuclear collisions in helium. In view of 
the theoretical interest in the stopping effect of 
still lighter or heavier gases, an examination of a 
large number of paired fragment tracks in 
hydrogen, deuterium, and xenon was performed 
and, at the same time, the former range measure- 
ments in helium and in argon were repeated in 
order to obtain more exact statistical data. 


EXPERIMENTAL METHOD 


The technical arrangement was the same as in 
the previous investigations, except for a few 
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Fic. 1. Paired fragment tracks in (a) hydrogen, (b) deuterium, (c) helium, (d) argon, (e) xenon. 


improvements. A 15-cm cloud chamber was used _ 


and, in order to reduce the stopping power of the 
vapor as much as possible, water was used as 
_ drop producing agent instead of alcohol-water 
mixtures. To avoid condensation of vapor on the 
uranium foils, 5 percent of NaCl was dissolved 
in the water gelatines. The stopping power of the 
gas mixtures was controlled by means of Po-a- 
particles from a source placed inside the chamber. 
The uranium layers were evaporated on }y-alu- 
minum foils and, in the case of xenon inves- 
tigations, on a mica foil. The thickness of the 
uranium layers was about 0.4 to 0.6 mg/cm*. For 
the calculation of the stopping power of the 
foils, 1.0 cm of air was taken equal to 1.5, 1.6, 
and 4.0 mg/cm? for mica, aluminum, and 
uranium, respectively. 

The tracks of fission fragments were thor- 
oughly examined and only those paired tracks 
were used for range determinations, where the 
end. points of both tracks were distinctly visible 
on the photographs. The end of the range, 
especially in light gases, is uniformly charac- 
terized by a tuft consisting of several short 
branches. Figure 1 demonstrates a few typical 
cases of paired fragment tracks in hydrogen, 
deuterium, helium, argon, and xenon. _Besides 
the typical characteristics of fission fragment 


tracks already described in previous papers* 
(heavy ionization, ejection of side branches, and 
fortuitous bending), slight differences appear in 
the extreme light and in heavy gases. As was to 
be expected from theory; the tracks in hydrogen 
and deuterium, as a rule, are rather straight from 
the start until pretty near the end of the range, 
whereas the tracks in argon and especially in 
xenon, due to numerous small collisions with 
heavy atoms, are submitted to a considerable 
bending over the greater part of the range. Also 
the length and the number of side branches are 
seen to change considerably with the gas in the 
chamber. In hydrogen, deuterium, and helium 
the tracks exhibit numerous branches; in argon, 
however, the branches usually are very short, 
and in xenon they are practically reduced to 
small lumps on the stem. Closer collisions which 
occur not infrequently in all gases and result in 
branches of considerable length are suitable for 


velocity determinations and will be discussed — 


below. 
COMPARISON OF RANGES 
The results of measurements of the reduced 
ranges* are given in the histograms of Fig. 2, 


4]. K. Béggild; K. J. Brostrém and T. Lauritsen, Kgl. 
Danske Vid. Sels. Math.-fys. Medd. 18, No. 4 (1940). 

* The reduced range rr is calculated according to the 
formula rr =38.0r/r_. mm, where r is the measured range of 


- 

wl 
dr 
of 
si 
du 
1 
th 
sol 
de 
in) 
vig shi 
sol 
= 
Xer 
Hel 
Hel: 
Hyc 
Deu 
4 — 
tr the 
cont 


gl. 


RANGE OF FISSION FRAGMENTS 283 
fy renge 

group 

2 (a) Xenon (ong range 


Number of tracks 


Fic. 2. Range distribution 
in different gases. 


where the groups of short and long range are 
drawn with opposite cross-hatching. The number 
of fission tracks in xenon and in helium is rather 
small, and the corresponding histograms repro- 
duce the numbers of tracks in range intervals of 
1 mm, whereas intervals of 0.5 mm are chosen in 
the other cases, where the statistical material is 
somewhat more copious. In hydrogen and 
deuterium, a few cases of very short range par- 
ticles were found. A comparison with other 
investigations suggests, however, that these very 
short range tracks originate in some accidental 
source of error as, for instance, the presence of 


TABLE I. Ranges in mm of air. 


Short range Long range Total 
Gas group group range 
Xenon 18 23 41 
n 
Rovtens data) 19 25 44 
n 19.4 23.9 43.3 
Helium 
(previous data) 23 30 53 
Helium 23 28 51 
Hydrogen 17.7 21.1 38.8 
Deuterium 18.9 22.5 41.4 


the tracks and rq the range of the a-particles from the 
control study. The constant 38.0 corresponds to the range 
of Po-a-particles at normal pressure and 15°C. 


water droplets; consequently, these tracks were 
not included in the statistics. The mean reduced 
range in mm of the two groups and the mean 
total range are given in Table I together with 
the. previous determinations of the mean ranges 
in helium and argon. 

Though the determination of the mean range 
is not too precise, the number of fission tracks 
being not very large, the values in Table | 
indicate characteristic differences in reduced 
ranges for the various stopping materials. If, in 
the first instance, we assume the electronic 
stopping power for fission fragments and a-par- 
ticles to depend in approximately the same 
manner on the atomic number of the gas, these 
differences must evidently be ascribed to the 
effect of nuclear collisions. The fact that nuclear 
collisions in the case of fission fragments, in 
contrast to the case of fast a-particles, actually 
have a considerable influence on the range is 
perhaps most clearly demonstrated by the 
results obtained in H and D. In these two gases, 
the electronic stopping must, of course, be 
identical, and the difference in the ranges gives a 
measure of the importance of nuclear stopping. 
It may be added that the longer range in D is 
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just what was to be expected, since, the H and 
D nuclei having the same charge, the lighter 
nuclei will receive more energy in correspond- 
ing collisions, except for the close collisions which 
are rare and unimportant for the mean stopping 
effect. 

In order to test in somewhat more detail the 
agreement with theory,? we note that the nuclear 
stopping power must be expected to be very 
nearly proportional to Z,?/M:2, Mzand Z2being the 
mass and charge of the nuclei of the atoms of the 
gas, and under the assumption that the de- 
pendence of the electronic stopping power on Z2 
may be taken from that of a-particles, we may 
therefore estimate the relative importance of the 
nuclear stopping in the various substances. This 
quantity should be directly connected with the 
reduced range. Thus, the reduced range is ex- 
pected to decrease with increasing ratio between 
nuclear and electronic stopping power. If, to a 
first approximation, we put the electronic stop- 
ping power proportional to Z,', the ratio of 
nuclear to electronic stopping, for heavier 
substances for which M;~2Zz2, is seen to increase 
roughly as Z;'. We should, therefore, expect the 
reduced range to decrease with increasing Z2 just 
corresponding to the experimental results in He, 
A, and Xe. 

It is obvious, however, that the marked short- 
ness of the reduced range in H and D compared 
with He is not compatible with the above as- 
sumptions concerning the electronic stopping 
power. In fact, even if on account of the large 
value of Z2/M; in H, these assumptions lead to 
a reduced range in H somewhat shorter than in 
He, it should still be longer than in Xe and, in 
particular, the reduced range in D should be 
expected to be even longer than in He. The 
measurements, therefore, seem to indicate a 
comparatively large electronic stopping power in 
H and D, pointing to an abnormally high effective 
charge of the fragments. It is of interest, how- 
ever, that this feature is just in accordance with 
theoretical considerations suggesting a much 
smaller probability of electron capture by the 
fragments in hydrogen than in other gases. 

A point which may require further discussion 
is the rather small difference between the reduced 
ranges in H and D, compared with the difference 
between the ranges in Xe and He. In fact, in D 


the relative importance of nuclear stopping is 
only half that in H, just as in He it is only about 
half that in Xe. A possible explanation would 
seem to be that, on account of the relatively 
large electronic stopping power, the nuclear 
stopping in H and D becomes of importance only 
for smaller fragment velocities than in other 
gases, in which case variations in the nuclear 
stopping would have a correspondingly smaller 
effect on the ranges. A more general discussion 


of these and other problems concerning the 


stopping of heavy nuclei will appear in a forth. 
coming paper by Bohr.® 

Extremely valuable information concerning 
the mechanism of uranium fission was given by 
Flammersfeld, Jensen, and Gentner,® and a com- 
parison with our range analysis might be of 
some interest. Unfortunately, however, the pos- 
sibility of a thorough analysis by means of the 
cloud-chamber method is rather limited in view 
of the paucity of statistical material, the finite 
thickness of the uranium layer, and a probably 


rather large range straggling. According to the _ 


theory,? the straggling of fission fragments is 
practically only due to the nuclear collisions, and 
the relative straggling is expected to be con- 
siderably larger in heavy than in light gases and 
smallest in hydrogen. Since the breadth of the 
range groups appears to be approximately the 
same in heavy and in light gases, it seems rather 
reasonable to attribute the main part of the 
spread in range to fluctuations in the initial 
velocity, charge, and mass of the fragrhents in 
the fission process. 

A comparison of the breadth of the range 
groups with the fluctuation in energy reported by 
Flammersfeld e¢ al. indicates that the range is 
less affected by variations in the fission process 
than is the energy of the fragments. This is just 
what should be expected from theoretical con- 
siderations. The appropriate formula given by 
Bohr? 

R=const. M,:Z;"!-V, 


where R is the range, M; and Z; the mass and 
nuclear charge of the fragment, and V the 
velocity, leads to the following equation, when 


5N. Bohr, Kgl. Danske Vid. Sels. Math.-fys. Medd. in 


press. 
*A. Flammersfeld, P. Jensen, and W. Gentner, Zeits. 
f. Physik 120, 450 (1943). 


FI 


| co! 

th 

has 

| ap) 

lig! 

atc 

| 

‘ 

is 


RANGE OF FISSION FRAGMENTS 285 


{long 


T T 
Hydrogen | 
‘ short range group 


| 
Deuterium 


{ short range group 
o + tong 


3s 


Fic. 3. Velocity-range relations 
in hydrogen and deuterium. 
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Z, is assumed to be proportional to Mi, 
R =const.- M,-"* 


where E is the energy. Since M,~’* is nearly 
constant for the fragments in the same group, 
the spread in range is roughly equal to half the 
spread in energy. 


VELOCITY-RANGE DETERMINATION FROM 
LARGE BRANCHES 


From quite general arguments regarding the 
dependence of stopping power on the charge, 
mass, and velocity of the incident particle, it 
appears evident that the heavy fragment group 
corresponds to the short range group, and the 
light fragment group to the long range group.’ 
Still, the verification of this assumption in 
previous experiments was based mainly on a 
statistical evaluation of the distribution of 
branches produced by close collisions with gas 
atoms,* whereas velocity determinations from 
individual large branchés only gave a mean veloc- 
ity range curve common to the two groups. 

As mentioned above, some of the fragment 
tracks in this work appear with individual large 
branches suitable for velocity determination and, 
consequently, it is possible to produce separate 
velocity-range curves for the two groups of 
fragments. Because of the fact that the fragment 

Fluegge), Kernphysikalische Tabellen (Verlagsbuchhand- 


lung. ulius Springer, Berlin, 1942). 
*N. Bohr, J. K. Boggild, K. ‘i Brostrgm, and T. 


Lauritsen, Phys. Rev. 58, 839 (1940 


6 & 16 
Range in mm. Range*in mm. 
tracks are comparatively straight in hydrogen 
and deuterium gas, the accuracy of the measure- 
ments of the angle between the direction of the 
branch is better here than in other gases. The 
velocity examination is, therefore, limited to 
close collisions with protons and deuterons. The 
velocity v of a fragment after a collision is given 
by the formula 


4 cos?0 


where k is the ratio of the masses of the gas 
nucleus and the fragment, u is the velocity of the 
gas nucleus, and @ the angle between the direc- 
tion of the stem before the collision and the 
direction of the branch. By collisions with 
protons and deuterons the ratio k is rather small 
and the following approximation can be used: 


v=u/2 cos @. 


The range of the branches was reduced to normal 
air conditions by means of the values for the 
stopping power of hydrogen relative to air, given 
by Gurney.® The velocity u for protons was ob- 
tained by means of the velocity-range relation 
for protons, given by Blackett and Lees,’® with 
the velocity increased by 5 percent," and for 
deuterons the velocity-range relation was found 


* R. W. Gurney, Proc. Roy. Soc. A107, 340 (1925). 
(93) M. S. Blackett and Lees, Proc. Roy. Soc. A134, 658 

ui M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, §95 (1937). 
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Fic. 4. Paired f ragment tracks in hydrogen. The up-going 
fragment track appears with two long proton branches, 
one close after the other, both suitable for velocity deter- 
mination. 


from the proton relation by a doubling of the 
abscissa. 

Figure 3 shows the result of ten individual 
velocity determinations in hydrogen and nine in 
deuterium. A rough estimation of the errors in 
the velocity determination due to the inaccuracy 
in the measurement of the angles and of the 
branch lengths is indicated by vertical lines 
through the points. The most probable initial 
velocities are calculated by means of the data 
given by Flammersfeld et al.; the values ob- 
tained are 9.010* and 13.3X10* cm/sec. for 
the heavy and the light fragment groups, re- 
spectively. Assuming the low initial velocity to 
correspond to the short range group, the ve- 
locity-range curves in both gases turn out to be 
rather reasonable and consistent with the 
previous results, whereas it seems to be quite 
impossible to produce reasonable relations if the 
low initial velocity were attributed to the long 
range group. Determinations of the velocity have 
not been carried out close to the end of the 
range, since the branches here are very short 
and, thus, such determinations would lead to 
considerable inaccuracy. Two very close col- 
lisions with protons, leading to extremely long 
branches, provide especially valuable velocity 
determinations rather early in the range of the 


short range group. The two branches in question 
are produced one close after the other by the 
same fragment, and a reproduction of the tracks 
is given in Fig. 4. 

Further argument to attribute the short range 
group to the heavy particle group is produced by 
a picture of the tracks of a close collision between 
a fission fragment and an argon nucleus, repro- 
duced in.Fig. 5. The down-going track has a 
range of 18.6 mm and the up-going track, in 
spite of the collision, a range of 22.1 mm. There. 
fore, the up-going track undoubtedly belongs to 
the long range group. The angles ¢ and 6* are 
found to be 22° and 45°, respectively, and the 
mass Mr of this fragment can be calculated by 
means of the formula 


Ma sing 


Mr sin(20+¢) 


Since the mass of the great majority of argon 
atoms is 40, the mass of the fragment in question 


Fic. 5. Paired fragment tracks in argon. The down-going 
fragment track is cut over by an almost parallel proton 
track, but the end point of the fragment tracks is clearly 
visible. The up-going long range fragment collides with an 
argon nucleus, and the mass of the fragment is found to 
be about 100 mass units. Hence, this long range fragment 
evidently belongs to the light fragment group. 


* » is the angle of deflection of the fragment and @ the 
angle between the direction of the argon branch and the 
direction of the stem before the collision. 
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is found to be 100+5 and, hence, this long range 
fragment evidently belongs to the light fragment 
group. 

Obviously, the above examination of indi- 
vidual large branches leads to velocity-range 
relations in conformity with previous, more 
indirectly produced information. 

Some recent experiments by Lassen™ have 
given, important information concerning the 
charge, energy, and energy loss of fission frag- 
ments passing through matter. In broad outline, 
Lassen’s experiments are in agreement with the 
cloud-chamber studies from this Institute; inti- 


2N. O. Lassen, Kgl. Danske Vid. Sels.: Math.-fys. 
Medd. 23, No. 2 (1945); Phys. Rev. 68, 142 (1945); Phys. 
Rev. 69, 137 (1946); Phys. Rev. 70, 577 (1946). 


mate conformity in all details is, however, not to 
be expected in view of the fact that close col- 
lisions take place chiefly in the later part of the 
range, whereas the accuracy of Lassen’s experi- 
ments is better in earlier parts of the range. 

The main part of the work described in the 
present paper was already carried out in 1943, but 
the completion and publication were delayed by 
wartime conditions. The authors wish to express 
their thanks to Professor Niels Bohr for his inter- 
est and advice regarding the experiments and 
their discussion. We are also indebted to the Nor- 
wegian Hydro-electrical Company which placed 
the xenon gas at our disposal and to Director S. 
Madsen, Copenhagen, for supplying us the argon 


gas. 


he 
ge 
en 
o- 
a 
In 
e- 
to i 
re = 
ne 
mn 
mn 


PHYSICAL REVIEW 


VOLUME 71, NUMBER 5 


Radioactive Isotopes of Re, Os, and Ir 


Leon J. GooDMAN AND M. L. Poot 
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(Received November 27, 1946) 


Two new radioactive isotopes have been found, Os'® produced by a Re(d,2n) reaction 
and Ir! produced by Ir(,2m) and by Os(d,m) reactions. The following characteristic values 


have been assigned to the isotopes in Re, Os, and Ir: 


Isotope HALF-LIFE 8-PARTICLE 7-Rays 
Reiss 92.8 hours 1.07 Mev No y-rays 
Reiss 18.9 hours 2.05 Mev 
Osi85 94.7 days 
Osi 31.9 hours 0.95 Mev 1.17 Mev 
10.7 days 0.091 Mev* 0.25 
Irisa 75 days 0.59 Mev 0.53 M 
Ir 19.0 hours 2.07 Mev 0.38 Mev and 1.65 Mev 


* Probably Internal Conversion Electrons 
The 17-day osmium activity has not been found. Chemical separations were carried out on 


the rhenium and osmium bombardments. 


INTRODUCTION 


N activity in rhenium was first found,’ when 
bombarding rhenium with neutrons from 
a Ra+Be source, to be of about 20 hours. With 
another Ra+Be source of neutrons there was 
found? a 20-hour and an 85-hour period. With 
fast neutrons® only an 18-hour period was found. 
The isotopes were placed and §-particle energy 
values were given after a deuteron bombardment 
of rhenium. The values assigned‘ were a 1.2-Mev 
B-particle for Re'** (90 hours) and a 2.5-Mev for 
Re!*8 (16-hour). The assignment of these isotopes 
was confirmed® by producing Re!®* by a W(d,2n) 
reaction. 

The radioactive isotopes in osmium have been 
reported as a 40-hour period‘ produced by slow 
neutrons and as 29.8-hour and 10-day periods® 
produced by neutron bombardment. Also re- 
ported’? have been a 32-hour period with a 
8-particle of 1.5 Mev and a 17-day period with a 
8-particle of 0.35 Mev. 

Two periods are known to be produced when 


M. Cork, and R. L. Thorton, es 
Rev. 52, 239 (19 7). 


«K. Sinna and H. Yamasaki, Phys. Rev. 55, 320 (1939). 
(1940) a H. Sullivan, Phys. Rev. 58, 276 


, Helv. Phys. Acta. 13, 219 (1940). 
Et aborg and G. Friedlander, Phys. Rev. 59, 400 


(1941). 


iridium is irradiated with neutrons, a 19-hour? 
period and a 60-day period.* The £-particle 
energy of the shorter period has been measured" 
and found to be 2.1 Mev. The periods were later 
reported" as 19.5-hour and a 68-day periods. 
These have been placed” using mass spectro- 


graphic technique at Ir™ and Ir’, respectively. - 


The investigation of the isotopes of rhenium, 
osmium and iridium was undertaken to clarify 
the existing disagreements and to give radio- 
active constants for the isotopes. 


Re'** AND Re!* 


Hilger Re metal (Lab. no. 9310) was activated 
by 10-Mev deuterons. ‘To the active sample was 
added 50 mg of Os metal carrier, which was then 
placed in an all-glass distilling apparatus. Con- 
centrated HNO; (30 cc) was added to dissolve 
the metals and the solution then boiled and 
brought nearly to dryness. More concentrated 
HNO; (20.cc) was then added and solution 
distilled until only a small amount of liquid (7 cc) 
remained in the distilling flask. The distillate, 
OsO,, was collected in alkaline (NH4)2Sz which 
precipitated the Os as OsSy The OsS, was 
filtered and measured for radioactivity. The Re 

8 E. Fermi, E. Amaldi, A. D’ tina, F. Rasetti, and 
E. Segré, Proc. Roy. Soc. A146, 483 (1934). 


inal Fomin and F. G. Houterman, Physik. Zeits. 


. Kamen, and S. Ruben, Phys. Rev. 
52, 375 “oh. 


4, R. Jaeckel, Zeits. f. Physik 110, 330 (1938). 
4 W. Rall, Phys. Rev. 70, 112 (1946). ; 


288 


Intenaity 


| MARCH 1, 1947 

| 

frac 
was 

rhe 

tha 

: frac 
T 

illus 

ee 1E. Amaldi, A. D'’Agostina, E. Fermi, B. Pontecorvo, be 

4 % F. Rasetti, and E. Segré, Proc. Roy. Soc. A149, 522 (1935). ae 

ne * 1. W. Kurtschatow, G. D. Latysctrew, L. M. Nemenow, 

tive 

The 

0.2 | 

shovy 

of 0. 

. 


RADIOACTIVE ISOTOPES OF RE, OS, AND IR 289 


Intensity 


i 


— 
= 
= 


Re+d 
Re fraction 
W not removed 


5 10 20 25 30 


= 
35 40 45 50 55 60 65 


Days After Bombardment 


Fic. 1. Decay of Re!** and Re'**, The intensities of the y-rays are shown as ten times their actual values. 


fraction was heated with a small amount (1.0 cc) 
of H.SO, until SO; fumes appeared. The solution 
was then cooled and diluted with H,O. The 
rhenium was precipitated as the sulfide, filtered 
and measured for radioactivity. It was supposed 
that any tungsten would remain in the: rhenium 
fraction. 

The decay curve of the rhenium fraction is 
illustrated in Fig. 1. The half-lives were found to 
be 18.9+0.2-hours and 92.8+0.2-hours. These 
values are assigned to Re'®* and Re!**, respec- 
tively, in accordance with previous workers.*® 
The large change in intensity over which the 
activity was followed made possible this claim of 
0.2 percent error. 

Aluminum absorptions of Re!** and Re'*®* are 


shown in.Fig. 2. The former gives a 8-end point’ 


of 0.442 g/cm? of aluminum or equivalent to 1.07 


Mev according to R=0.526E—0.094." 
gives a B-end point of 0.995 g/cm? of aluminum, 
equivalent to 2.05 Mev. There are no y-rays 
associated with Re'®* and the number of y-rays 
associated with Re'** are less than one per four 
B-particles. 


The relative saturation intensities are 6800: 


18000 for Re'**: Re'**. The previously reported | 


Re(d,a) reaction producing W'® was expected 
but not found. Assuming either W'® or Re'* 
(produced by stray neutrons) were present in 
the tail end of Fig. 1, the relative saturation 
intensity would be the order of unity. 


Os!*! 


In order to investigate the decay of osmium, a 
0.5 gram sample of osmium metal was irradiated 


4 B. W. Sargent, Can. J. Research 17A, 82 (1939). 
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Fic. 2. Absorption in aluminum of §-rays from, Re™* 
and Re!#, 


with slow neutrons produced by Be+ (5 Mev) 
for 17 hours. Figure 3 shows the decay of this 
sample. 

The 31.9-hour period is seen as the principal 
activity present. The previously reported 17-day 
period has not been found even though long 
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Fic. 3. Decay of Os’ produced by slow neutrons. 
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Fic. 4. Absorption in aluminum of 8-rays from Os, 


bombardments were made. Assuming the 17-day 
period present in the tail end of Fig. 3 the ratio 


of the saturation intensities found would be of 


the order of 250: 1. 

Figure 4 shows an aluminum absorption during 
the 31.9-hour period of the 8-rays from an os- 
mium fraction of an Os+d bombardment. The 
B-end point of 0.400 g/cm? of aluminum, equiva- 
lent to 0.95 Mev, was found. A lead absorption 
revealed a y-ray of half-thickness of 0.40 inch or 
of 1.17 Mev. 

A weak activity of about 30-hour half-life has 
been obtained by bombarding iridium with 
deuterons, presumably caused by an Ir(d,a) 
reaction leading to Os'. 


Og!85 


A new period has been found in the osmium 
fraction after a deuteron bombardment of Re. 
The chemical separations were carried out as 
previously described. The decay of this isotope 
is shown in Fig. 5 and is seen to have a half-life 
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Transmutations in The Os Region 
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of 94.7+2.0 days. The isotope is predominately 
y-ray emitting and has been tentatively placed 
at Os'®. It is believed that Os'® decays by 
K-electron capture. The previously reported" 
30-hour isomer of Os!®, produced by Re(d,m) was 
not confirmed when very pure rhenium was used. 
In Fig. 6 is shown a chart which summarizes 
the reactions and values given in this paper. 


A new period has been found in iridium. It 
has been made by irradiating iridium with fast 
neutrons and has been found in the iridium 
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= 20 
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chemical fraction after bombarding osmium with 
deuterons. The decay of the iridium activity 
induced by neutrens from Li+d is shown in 
Fig. 7. The half-lives are seen to be approxi- 
mately 10.7 days and 75 days. 

Since the 10.7-day period has been made by 
fast neutrons on iridium and by Os(d,m) and has 
not been found with slow neutrons or with 
deuterons on iridium, the activity is therefore 
ascribed to Ir!®. From the various bombard- 
ments the half-life of 10.7++0.3 days is associated 
with this period. 

Associated with Ir! is an internal conversion 
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Fic. 7. Decay of Ir! and Ir. 


*L. J. Goodman and M. L. Pool, Phys. Rev. 70, 112 (1946). 


{ 
day 
atio 
ring 
he 
The 
iva- | 
tion 
h or 
has || 
vith 
d,a) | 
5000 
ium 2000 
Re. 
t as 
-life 
| 
| 


292 L. J. GOODMAN 
100 
a 
2 
4 
z 
2.07 Mev due to Ie™* 
os H 


20 2 Seoie 


Fic. 8. Absorption in aluminum of 8-rays of Ir'* 
and Ir™. 


electron or very soft 8-particle and a y-ray. The 
electron has an end point of 0.0110 g/cm? of 
aluminum or an energy of 91 Kev. The y-ray 
when absorbed in lead, has a half-thickness of 
0.039 inch or an energy of 0.25 Mev. X-rays have 
also been found to be associated with this 
period, suggesting K-electron capture. 


AND Ir! 


The 19.5-hour and 68-day period" in iridium 
have been placed” at Ir! and Ir!, respectively. 
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Half-life measurements of these isotopes have 
given values of 19.0+0.2-hours for Ir!™ and of 
75+3-days for Ir’ agreeing with the Previously 
reported periods. 

The 8-end point of Ir', as shown in Fig, 8 
has been found to be 0.216+0.007 g/cm? of 
aluminum, equivalent to 0.59+0.01 Mev. The 
Ir'™ 8-end point of 1.00+0.02 g/cm? of alumi. 
num, equivalent to 2.07+0.03 Mev, which 
agrees well with the previously reported" valye 
of 0.970+0.030 g/cm’. 

The y-rays of Ir’ and Ir’ were previously 
reported’*® as 1.35 Mev and 0.63 Mey, respec. 
tively. The values of 1.65 Mev and 0.38 Mev for 
Ir'™ and 0.52 Mev for Ir'® have been obtained, 
. The procedure devised to study the y-radi- 
ation was as follows. A number of lead and tin 
absorption curves were. made at intervals after 
an Ir+d bombardment and the intensities of the 
y-rays, when extrapolated back to zero thickness 
of lead (Jo), weré plotted as ordinates against 
time after bombardment. An absorption curve 
illustrating how these values of J) were obtained 
is shown in Fig. 9. This absorption curve shows 
a y-ray of half-thickness in lead of 0.500 inch, 
equivalent to 1.65 Mev, and a y-ray of half- 
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Fic. 9. Absorption in lead of y-rays from Ir+d 
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18 C, M. Witcher, Phys. Rev. 60, 32 (1941). 


16C. E. Mandeville and H. W. Fulbright, Phys. Rev. 
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64, 265 (1943). 


Oo we 


paar ay 


ig 

a: 

} 

{ 

| 

4 

‘ 

Te — 

\ ~ 

\ 

= \ 

\ ~.390 Mev 

< 

| 


RADIOACTIVE ISOTOPES OF RE, OS, AND IR 293 


Decay Of 7-Rays 
from Ir +a Bomb 
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Fic. 10. Decay curve of Io. The values of the half-thicknesses in tin and lead are shown for 
each absorption. 


thickness in lead of 0.105 inch, equivalent to 
0.39 Mev. The decrease in intensity of these zero 
thickness values with time is shown in Fig. 10. 
The J» of the y-ray of 1.65 Mev is seen to decay 
on an 18.8-hour half-life. 

The decay curve of the other J» is seen to have 
two half-life components, a 19-hour period and a 
75-day period. It was found that the half-thick- 
ness value in lead of the curve associated with 
this J) increased with each absorption after the 
first, varying in half-thickness value from 0.105 
inch for the first absorption curve to a value of 
0.163 inch when the 19-hour component of the 
decay curve was no longer present. In order to 
explain the fact that there are different energies 
for each absorption, it is necessary to assume 
that there are two y-rays, one associated with 
each period, which cannot be resolved by lead 
absorption technique. Using the half-thickness 
in lead value of 0.163 inch equivalent to 0.52 
Mev, as that of the energy of the y-ray associated 
with Ir’ (75-day), it has been calculated that 
the half-thickness of the y-ray associated with 


Ir! (19-hour) is 0.102 inch of lead, equivalent 
to 0.38 Mev. 


POSSIBLE ENERGY LEVELS FOR Ir'** 


Ir has been found to have a §-particle of 
2.07 Mev and y-rays of 1.65 Mev and 0.38 Mev. 
The two y-rays are of approximately the same 
intensity and there are in the neighborhood of 
four 8-particles for each group of the two 7-rays. 
This suggests that there may be two alternate 
methods of decay. First by simple emission of a 
2.07-Mev §-particle, and second by consecutive 
emission of the two y-rays, and then by emitting 
a low energy §-particle (about 0.03 Mev). This 
low energy §-particle has not as yet been found. 
On the basis of the above proposed scheme of 
decay the first procedure takes place 80 percent 
of the time. 

The authors are grateful for the support re- 
ceived from Mr. Julius F. Stone, The Ohio State 
University Development Fund and The Ohio 
State University Research Foundation. 
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Theory of the Refraction and the Diffraction of Neutrons by Crystals 
M. L. GotppercEr, Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
AND 


FREDERICK SE1Tz, Physics Department, Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
(Received November 27, 1946) 


The equations for the elastic scattering of neutrons by a single crystal, governing simple 
refraction, reflection, and Laue-Bragg scattering, are derived under the assumption that the 
scattering and absorption cross sections are independent of spin. A brief summary of the 
results is given in the final section of the paper. 


1. INTRODUCTION 


HE experiments on the scattering of neu- 
trons by single crystals carried out at the 
Argonne Laboratories by Fermi and Zinn and at 
the Clinton Laboratories by Borst and Wollan 
have given conclusive evidence that an appre- 
ciable part of the scattered intensity is coherent. 
In fact, a substantial part of the scattering takes 
place in accordance with the Laue-Bragg equa- 
tions. The purpose of the present paper is to 
investigate the theoretical background for the 
scattering somewhat more fully than has been 
done previously! and to summarize the results in 
a form that may be of use in the course of the 
development of the experimental work. 

The scattering of slow neutrons by crystals has 
been of considerable interest since the earliest 
days of neutron physics. Following Fermi’s 
development? of a simplified or semi-empirical 
method of treating the interaction between 
neutron and nucleus, Wick*® showed that one 
should expect highly crystalline media to scatter 
slow neutrons in a way that is very different 
from that expected for gases or completely 


1 This document is based on work performed under Con- 
tract No. W-35-058-eng-71 for the Manhattan Project. 
Part of the information covered in this document appeared 
in Report CP-2419, and a more —e. survey will appear 
in volume IIIB, Division IV of the Manhattan Project 
Technical Series, as part of the contribution of the Clinton 
Laboratories. Some of the results derived in this report, 
particularly those for the index of refraction, were derived 
earlier by Fermi for simple cases. Fermi’s measurements of 
the total reflection of thermal neutrons by graphite and 
subsequent measurements of Bragg scattering by Zinn 
and Borst furnished the incentive for much of the work 
described here. Professor W. E. Lamb has informed us 
that he investigated theoretical aspects of the problem of 
neutron refraction as early as 1940. 

2 E. Fermi, Ricerca Scientifica, 7, Part 2, 13 (1936). 

*G. C. Wick, Phys. Zeits. 38, 403 (1937). 


amorphous materials because of the interference 
effects which occur. Wick’s work has formed the 
basis for much of the theoretical work in this 
field. Following a very similar line of reasoning, 
Teller* pointed out that one should expect 
ortho- and parahydrogen to possess markedly 
different scattering cross sections because of 
interference, provided the large cross section for 
the scattering of slow- neutrons by protons could 
be ascribed to a virtual singlet level of the type 
used in the Breit-Wigner formalism of resonance 
scattering. The subsequent experimental veri- 
fication® of the quantitative predictions of Teller 
and Schwinger‘ on the basis of Teller’s original 
suggestion demonstrated that one can expect to 
obtain a reasonable description of the scattering 
of slow neutrons by polyatomic systems with the 
use of wave mechanics and the Breit-Wigner 
formalism. 

In the period following Wick’s work, numerous 
investigators extended his treatment of the scat- 
tering of neutrons by crystals. Most prominent 
among these are the work of Pomerantschuk,* 
Van Vleck,’ Halpern, Hamermesh and Johnson, 
Seeger and Teller,® and Weinstock.'° Pomerant- 
schuk examined more carefully than Wick the 
influence of low temperatures upon the scattering 
cross section. Van Vleck investigated the scatter- 


4 J. Schwinger and E. Teller, Phys. Rev. 51, 775 (1937); 
J. Schwinger and E. Teller, Phys. Rev. 52, 286 (1937). 

5 J. Halpern, I. Estermann, and O. Stern, Phys. Rev. 
52, 142 (1937); L. W. Alvarez and K. S. Pitzer, Phys. Rev. 
55, 596 (1939). 

1938) Pomerantschuk, Phys. Zeits. Sowjetunion 13, 65 

7 J. H. Van Vleck, Phys. Rev. 55, 924 (1939). 

80. Halpern, M. Hamermesh, and M. H. Johnson, 
Phys. Rev. 59, 981 (1941). 

*R. J. Seeger and E. Teller, Phys. Rev. 62, 37 (1942). 

10R. Weinstock, Phys. Rev. 65, 1 (1944). 
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ing of neutrons by paramagnetic media. Halpern, 


Hamermesh, and Johnson considered the elastic . 


or Laue-Bragg scattering of neutrons by a simple 
lattice. Seeger and Teller treated in a qualitative 
fashion the consequences of the selection rules 
which govern inelastic collisions of neutrons with 
a lattice in the approximation in which one 
thermal quantum is emitted or absorbed in the 
process. Finally, Weinstock presented a more 
complete treatment of the theory of inelastic 
collisions and applied his results to a computation 
of the influence of temperature on the scattering 
cross section of a polycrystalline medium as a 
function of temperature. His work was carried 
out in the approximation in which each crystal 
scatters a negligible portion of the intensity of a 
monochromatic beam of neutrons which im- 
pinges upon it. Use of the foregoing theoretical 
work will be made in the following presentation. 

It should be emphasized that a number of 
experiments which were carried out prior to 
those at Argonne and Clinton gave strong indica- 
tions that polycrystalline media scatter differ- 
ently from non-crystalline ones. In this con- 
nection should be mentioned the work of Beyer 
and Whitaker" on the scattering of neutrons by 
salts, that of Anderson” and co-workers on the 
penetration of very slow neutrons through thick 
layers of graphite, and the experiments of Rain- 
water and Havens" on the variation with neutron 
velocity of the scattering cross section of solids. 

The present paper will treat the theory of 
refraction and diffraction from several points of 
view and will derive equations both for the 
complex index of refraction and for the reflec- 
tivity when the Laue-Bragg conditions are 


satisfied. For simplicity it will be assumed that | 


the scattering cross sections are independent of 
spin orientation. This condition is satisfied when 
the scattering nuclei have zero spin, as is true 
for the isotopes of even atomic mass above N". 
There is considerable evidence that the equations 
based on this assumption are reasonably accurate 
for many nuclei having a finite spin. For ex- 


1H. G. Beyer and W. D. Whitaker, Phys. Rev. 57, 976 
oom. See also the more recent work by F. C. Nix and 

. F. Clement, Phys. Rev. 68, 159 (1945). 

® Anderson, Fermi, and Wood, CP-718, CP-781 (Man- 
hattan District Reports). 

and Havens, CP-1962 (Manhattan District 


ample, the alkali halides exhibit good Bragg 
reflection. 


2. EQUATIONS FOR A SCATTERING MEDIUM 


According to conventional collision theory,“ 
the wave function describing an incident plane 
wave of particles and the waves scattered by a 
spherically symmetric center of force is 


f(6)/r, (1) 


at large distances from the center, where k is the 
wave number vector of the incident beam of 
particles; r is the positional vector; and 


4 


The center of force is assumed to be at the origin 
of coordinates. In Eq. (2), P, is a Legendre 
polynomial of order ; n(m) is the shift in phase of 
the radial wave function associated with P, 
relative to the phase in the absence of the force 
field. This phase shift is real if the center does 
not absorb; however, it is complex in the more 
general case in which the center may absorb the 
particles as well as scatter them.'® In order to 
treat the case in which (m) is complex we shall 
write it in the form 


n(n) = ne(m) +in(n) (3) 


in which both 9,(m) and are real. 

An analysis of the amplitude of the scattered 
wave in (1) shows that the scattering cross 
section o, of the center is given by 


—29i(n 
n(2n+1)e 


(4) 


Similarly an examination of the net flux of 
current through a sphere surrounding the center 
shows that the absorption cross section is 


2 
sinh2n(n). (5) 


4 See for example N. F. Mott and H. S. W. Massey, The 
He of Atomic Collisions (Oxford University Press, 

4b We are indebted to Professor Schwinger for an in- 
formative discussion of complex phase shifts. 
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The total cross section, ¢,=0,+0¢, is given by 
the relation 


(6) 


in which Im implies “imaginary part of.”” In the 
present case 


2 
(7) 


In ordinary problems of nuclear scattering the 
phase shifts are small, corresponding to the fact 
that cross sections are usually small compared to 
a/k? (about 2.6 cm? or 2.610? barns for 
thermal neutrons). In this approximation Eqs. 
(4), (5), and (7) become: 


4 

nr 

+1) +940): (9) 


+ —nr(n)). (10) 


We readily see from these equations that 7;(0) 
is very small compared with 7,(0), except when 
o, is much larger than o,. The ratio of the two 
components is 


The components of phase become comparable 
only when go, is of the order of 10‘ barns for 
ordinary values of 

If we expand the wave function (1) in the 
vicinity of r=0, we find that in the approxima- 
tion in which f(@) =/k the leading terms are 


Y=n/kr+1+---. (11) 
Thus near the origin y in general has the form 


(12) 
where 
ka_1/ao=. (13) 
In view of Eqs. (8)—(10), this means that 
ka_; koa 
ao 4n 
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3. WAVE FUNCTIONS IN A CRYSTAL 


We are now in a position to examine neutron 
wave functions in a crystal. We shall attempt to 
find a function of the form 


=x(r)e*'*, (15) 


in which x has the periodicity of the lattice, r is 
the positional vector, and k’, is the wave number 
vector. For simplicity we may assume that the 
crystal is a simple cubic lattice although the 
results which we shall obtain are easy to gener- 
alize to a more complicated case. We shall also 
assume that the energy E of the wave function 
is related to the wave number k outside the 
crystal through the relation 


E=h*k?/2m. (16) 


Since the wave function possesses a 1/r singu- 
larity of the type (12) at the position of each 
nucleus in the lattice, and since the potential of 
the neutron can be taken to be zero in the region 
of space between nuclei, it follows that the equa- 
tion satisfied by the wave functions within the 
lattice is 


Ay+k*y = —4ra_, >, (17) 


The right-hand side of this equation is a sum of 
delta-functions over the positions r, of the 
various nuclei, appropriately modified by a phase 


. factor in accordance with the form of (15), 


because the Laplacian applied to 1/r gives 
—4ri(r). The corresponding form of the equation 
for the function x(r) is 


V2x + - gradx — (k’?—k?)x 
= —4ra_; 6(r—f,). 


Since x is a periodic function, it can be expanded 
in a Fourier series 


(18) 


x= Dx (19) 
The summation is to be carried out over all 
vectors K of the reciprocal lattice, that is, over 
all vectors the components of which are integer 
multiples of 2x/a, where a is the lattice constant. 


16 The mathematical approach presented in this section 
was proposed to us by eer ea The results ob- 


tained with this method were first obtained by the some- — 


what less rigorous method described in Section 5. 
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If we introduce (19) into (18) and expand the 
sum of the delta-functions into a Fourier series, 


we obtain 
cx( —K?— 2k’-K —k’?+k’) —4nra_;/v, (20) 


where v=a? is the volume of the unit cell of the 
lattice. Thus when we solve for cx and substitute 
into (19) and (15) 


4ra_; exp[i(k’+K)-r] 


The series converges very slowly ; in fact its con- 
vergence behavior is the same as for the series 
expansion of the function of 1/r. It is convenient 
for subsequent use to place (21) in the form 


i(k’ . 
exp[i(k’ +) r] 


dr(y). 
k’ 2_ 
(k’+y) (22) 


It now remains to relate a_; to the physically 
interesting quantities o, and a, (see Eq. (14)). In 
order to use Eq. (14) we must first determine do 
from (22). To do this, we subtract from (22) the 
‘quantity 
ras 

r (27) *y? 
exp[i(k’+y) -r] 
=4 f 
(23)*(k’+y)? 


dr(y) (23) 


and obtain 


r 


— f dr(y) expLi(k’ 


é(u—K 
(u—K) (24) 


This expression converges for r=0, and its 
limiting value is 


5(y—K) 
e= dr 
f Obs 


(25) 
(20)*(k’ +p)? 


If (14) is used to replace the ratio a_;/ao, this 
equation may be employed to determine the 


relation between k’, k, and the observed cross 
sections of the nucleus involved. 

In the remainder of this section we shall 
restrict our attention to the case in which k” and 
k? are very nearly equal, as will be true whenever 
the phase shifts are small, and in which (k’+K)? 
and k* are nearly equal only when K=0. This is 
the case in which the Laue-Bragg conditions for 
diffraction are not satisfied so that the crystal 
behaves like a simple optical medium. We shall 
consider Bragg reflection in detail in the next 
section. Under these conditions the first term in 
the sum in (25) is much larger than the other 
terms appearing behind the integral sign. As a 
result the equation can be written 


ao 4r 1 


ay wv k’?—k?* 
or 


in which mo is the number of atoms per unit 
volume. 

It is interesting to note that k’ will be complex 
in the general case in which ¢, does not vanish. 
Thus in this case the functions (15) have the 
form of attenuated waves. 

If o, is zero, Eq. (27) becomes 


k’?—k? = (28) 


(26) 


The index of refraction m of the crystalline 
medium, defined by the relation 


k’=nk (29) 
is given by the equation 
k? 8x*mE 


It is to be noted that there is an arbitrariness 


in sign before the square root in Eqs. (14), (28), 


and (30). This sign cannot be determined merely 
from knowledge of the scattering cross sections, 
but requires additional information. According 
to the Breit-Wigner formalism of resonance 
levels, the sign which appears before the root 
depends upon the distribution of resonance levels 
of the compound nucleus relative to the zero 
energy represented by the energy of the incident 
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neutron plus that of the nucleus. Resonance 
levels which have higher energy contribute to the 
positive sign, those with lower energies to the 
negative sign, just as in the corresponding optical 
case. 

In cases in which the nuclear scattering cross 
section is of the order of 10-*4 cm?, as is commonly 
the case, n?—1 is of the order of 10~* for thermal 
neutrons (k~5X10* Thus the index of 
refraction differs from unity by a very small 
amount, in close analogy with the situation’ for 
x-rays. 

The wave functions of the form (15) may also 
be derived with the use of the cellular approxima- 
tion, which has proved very useful in the treat- 
ment of electron waves in solids, in the case in 
which the phase shift is real (i.e., o,=0). A 
straightforward application of this procedure 
shows that as long as the phase shift is small, the 
form of x within any given cell of the lattice is 
accurately’ represented by the function 


xo=1+10/kor (31) 


in which ko is the value of k for which k’ is zero, 
that is, according to (28), 


ko? = +no(410,) 
and 4 is the phase shift for k=ko, that is, 
no= + 


This form of x is valid for values of k’ outside 
the first zone, provided k’ is not too close to the 
surfaces of the Brillouin zones. The same type of 
calculation shows that the effective neutron 
mass inside the crystal is the actual mass. This 
conclusion is, of course, apparent from Eq. (27) 
in the approximation in which o, vanishes since 
the coefficients of k? and k” are equal. 

Equation (31) may also be represented in the 
form 


in which is the “extrapolation distance’ of the 
wave function near the origin, which is related 


17 In the in which the cells of the crystal 
lattice are replaced by spheres, as is possible in reasonably 
ones lattices, the wave function xo will have the 
orm 

xo= (1/kor) sinko(r —10). 


xo=1—10/r, (32) 


to the other quantities by the equations 
Roto=—n0; ko? =4amoro; (33) 


Since (¢,/4)* is of the order of 10-* for most 
nuclei, it follows from the last of the relations 
(33) that the second term on the right-hand side 
of (32) is negligible when r corresponds to dis. 
tances midway between atoms. In other words 
the function x» is almost constant everywhere 
except in the immediate vicinity of the nucleus 
when the scattering cross section is of the order of 
10-*4 cm?. 

When 1» is used in place of (¢,)', the equation 
for the refractive index n is 


n?—1= —ky?/k?= —Annory/k*. (34) 


In the event that there is more than one atom 
per unit cell, the more general form of Eq. (17) is 


= > exp(ik’-r)a_i(j) 
X6(r—1rj;—1,). (35) 


Here the vector r, is summed over the cells of the 
lattice whereas r; represents the positional coor. 
dinate of the jth atom in the unit cell and is 
summed over the atoms in this cell. a_:(j) is the 
coefficient of 1/|r—r;| in the expansion of the 
wave function about the point r; We shall let 
ao(j) designate the corresponding value of the 
constant term in the expansion about r;. The 
relation between do(j) and a_:(j) analogous to 
(14) is 


ka_1(j) Raa! 
ao(J) 


=ké(j), (36) 


in which o,/ and o,/ are the scattering and ab- 

sorption cross sections for the jth nucleus. ° 
The equation for cx in the expansion of (19), 

analogous to (20) is readily found to be 


cx[k?—(k’+K)?] 
= —Anny a_(j) exp(—iK-r)), (37) 
in which mo is the number of unit cells per unit 


volume of the crystal. Similarly the family of 
equations, analogous to (25), relating the ao(j) 


is 
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and a_;(j) are 


(38) 


f 


The vector p has the same meaning as in (25). 

As previously, only the first term in the sum 
(38) is large when k” and k* are very nearly 
equal (and (k’+K)? is not nearly equal to k? for 
any K different from zero). Thus in this case Eq. 
(38) may be written 


ao(t) >; a_-1(j) /(R’?—k?). (39) 


Since the right-hand side of this is independent of 
i, it follows that 


do(t) =C, (40) 


in which C is a constant independent of 7. Using 
this relation, we find that (39) may be placed in 
the form 


which replaces Eq. (27). The new form of y is 


exp(tk’-r) 


(41a) 


Equation (41) may also be derived with the 
use of the cellular approximation in the case in 
which the a,‘ are zero. In this case one surrounds 
each nucleus with a polyhedral cell and deter- 
mines the wave functions within the cell by 
means of the usual boundary conditions. As 
previously, it is found that the wave function has 
the form (32) in each sub-cell of the unit cell of 
the lattice; that is, the function x is constant 
everywhere except in the immediate vicinity of 
the nuclei. 

The fact that x is constant at regions not too 
close to the immediate vicinity of the nuclei has 
very important consequences when the solid is 
treated like an optical medium. In the conven- 
tional treatment of the refracting properties of a 
continuous optical medium, the relations between 
incident, transmitted, and reflected rays are 
determined by fitting boundary conditions at the 
surface of the solid, which in the simplest cases 
is taken to be a plane. Since x is constant at 


Xx, ; 6(y—K)a_.(j) ] a_,(t) 
(k’+y)?—k? 


no(2m)*(k’ 


distances of the order of 10-* cm from the nuclei 
in the present problem, it follows that the neutron 
wave functions may be taken to have the form 
exp(ik’-r) at a mathematical surface of the solid 
which lies at the periphery of the outermost layer 
of atoms of the crystal. If this surface is reason- 
ably flat in the sense that large areas of the 
crystal surface are composed of the same con- 
tinuous crystallographic plane, the conventional 
results relating index of refraction, reflection, and 
transmission of a continuous optical medium, 
derived for the planar case, can be employed. 
For example, in the case in which the index of 
refraction is real, the reflectivity is given by the 
relation 


(n?—sin26)!— =] 
(n?— sin?6)'+-cosé 


in which @ is the angle of incidence of the incident 
beam, that is, the angle made with a normal to 
the surface. If m is greater than unity this remains 
less than 1; however, if m is less than unity, the 
reflectivity is unity for values of @ greater than 
those for which 


sin@>n. 


Since m is very near to unity in the cases in which 
we are principally interested here, we may write 
n=1—p and use the. relation 
cosé =singy. The reflectivity then becomes 


1—[1—(2u/¢*)]*)? 


Perfect reflectivity occurs for values of ¢ less 
than (2u)*. For typical nuclear cross sections, 
one finds that the angle at which the reflectivity 
becomes total is of the order of 10 minutes of arc. 
This type of reflectivity was first sought for and 
found by Fermi and Zinn. 

As in the conventional optical case, the com- 
ponent of k’ in the direction normal to the 
crystal surface and inside the solid becomes 
imaginary when sin@>n, if m is less than unity. 
That is, in the region of total reflection the 


(42) 
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transmitted wave is dampened exponentially. 
The depth of penetration d for which the am- 
plitude falls to 1/e times its value at the surface is 


i 
k(n?—sin?0)* k(2u)* 


in which \ is the wave-length. We find that d is 
of the order of 10°A for typical values of yu, that 
is, of the order of 100A. This distance is suf- 
ficiently small that experimental determinations 
of the existence of total reflection will be seriously 
influenced by the presence on the surface of a 
contaminating layer of the order of 100 atom 
distances in thickness. For example, the oxide 
layer on most metals will have an important 
influence on the reflectivity. 

The recent experiments of Fermi,!* Zinn, 
Marshall, Borst,” and Wollan, and the somewhat 
earlier work of Goldberger carried out in coope- 
ration with ‘Zinn, seem to show that most nuclei 
scatter with negative phase shift, corresponding 
to index of refraction less than unity. (Man- 
ganese appears to be the only exception dis- 
covered to date.) This result seems to imply that 
for most nuclei the influence of stable neutron 
levels having negative energy on the dispersion 
properties of nuclei for thermal neutrons is 
greater than that of the virtual levels which lie 


above the thermal region. If further work jus- ~ 


tifies this conclusion we need not be too sur- 
prised: Practically all transparent media have 
an index of refraction greater than unity in the 
visible portion of the optical spectrum ; it would 
hardly seem more startling to find a similar 
uniform principle operating for thermal neutrons 
in cases in which the scattering nucleus does not 
possess a virtual level lying close to the thermal 


region. 
4. LAUE-BRAGG REFLECTION 


We are now in a position to discuss the 
analog of ordinary x-ray diffraction. From the 
mathematical viewpoint, this type of diffraction 
occurs when k’? and (k’+K)* become equal for a 
value of K different from zero. We shall designate 


18 E, Fermi, W. H. Zinn, and L. Marshall, Abstracts E2, 
E4, and E5, Chicago Meeting, American Physical Society, 
June, 1946. 

19 As yet unpublished work. 


this value of K by K’ for the purposes of the fol- 
lowing discussion. When this condition is satis. 
fied, two terms in the sum under the integral sign 
in Eq. (38) are large instead of only one, namely 
that for which K is equal to 0 and to K’. Thus 
the equation relating k and k’ which replaces 
(39) is 


a-1(J) 
k’2?—k? 
@-1(j) exp(—iK’-r,) 


+exp(iK’-r,) ). 


ao(t) 


If we write 


(k’ +K’)?—k? 


Eq. (43) may be placed in the form 
ao(t) = (46) 


When the value of a_:(j) derived from (46) is 
substituted into (44) and (45), the following two 
simultaneous equations result: 


+C2 exp GK’ 
Ci 
D; 


The determinantal compatibility equation asso- 
ciated with (47) gives the relation between k* 
and k’*. We may expect the roots of this equation 
to exhibit all of the characteristic features of the 
theory of Brillouin zones in the case in which the 
£(z) are real. That is k? will be single valued and 
real for all real values of k’ except at the zone 
surfaces for which k” and (k’+K’)? are equal. 
Here k*® will be real and double-valued. Con- 
versely, in the neighborhood of the zone surface 
there will be two real values of k” associated 
with each value of k?. These values of k” will 
be complex when the £(z) are complex. In other 
words, there are two sets of functions y associated 
with each value of k? when the Laue-Bragg con- 
dition is satisfied. From (19), (37), and (46) we 
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find that these functions have the form 
C, exp[sK- (r—r,)] 
(k’ + K)?—k? 


exp(tk’- 


+ C: exp(tk’ - r) 
K 


in which C,; and C, are the solutions of Eqs. (47). 
We see that these functions may be regarded as 
the superposition of two waves traveling in the 
direction k’ and k’+K’, respectively. The rela- 
tive amplitude with which each solution appears 
in the final form of the wave function within the 
crystal is to be determined by the boundary con- 
ditions. These in turn depend upon the physical 
arrangement in a given experiment. For example, 
if a neutron beam is diffracted in such a way that 
it enters and emerges through the same surface 


(the case of Bragg reflection), the solution within’ 


the crystal must join the incident and emergent 
solutions at the surface of incidence. In addition, 
it must satisfy the condition that there be no 
incident wave on the rear face of the specimen. 

Fortunately, we may note that Eqs. (47) 
are identical with those which occur in the 


‘dynamical theory of x-ray diffraction” if C, and 


C, are interpreted as being equivalent to the 
amplitudes Do” and Dx’, respectively, of the 
displacement vectors which are polarized in a 
direction normal to the plane containing k and 
k’, and if the quantities 


mo &(j) and mo &(j) 


occurring as coefficients of C, and C; in the right- 
hand side of (47) are regarded as being the 
equivalent of k® times the Fourier components 
of the complex polarizability of the medium. 
Moreover, it is conventional to have Do and Dx: 
satisfy exactly the same boundary conditions as 
the corresponding amplitude constants in our 
problem do. Thus we may use all of the results 
of the dynamical theory of x-ray diffraction for 
the case of the normally polarized wave, after 


*W. H. Zachariasen, X-Ray Diffraction in Crystals 
(Wiley and Sons, New York). See . 9, Chapter 3 in 
perSetee. Weare deeply indebted to Professor Zachariasen 
or numerous valuable discussions of the material in this 
paper and for the use of an advanced copy of the manu- 
script of his book. 


(k’+K)*—k’ 


| (48) 


suitably readjusting the constants which appear 
in the equations. 

It is readily seen that Eq. (48) may be placed 
in the form 


¥ exp(tk’ -r)xx-(1) 
+C, (48a) 


in which xy and xx-4x’ are identical functions of 
the variables k’ and k’+K’, respectively. x’ in 
(48a) is, in turn, defined by the same series as the 
corresponding function xx, appearing in (41a), 
which we have seen to be nearly constant every- 
where except in the immediate vicinity of the 
nuclei. It follows that as long as r is not in the 
immediate vicinity of the nuclei (48a) may be 
written in the approximate form 


=A, exp(ik’ -r)+ A: exp[i(k’+K’)-r], (48b) 


where Ci/C2. 

Now in satisfying conditions of continuity for 
incoming and transmitted waves in the dynamical 
theory of x-ray diffraction, it is commonly 
assumed that the wave inside the crystal has the 
form (48b). This assumption is highly accurate 
in the present case near the mathematical 
boundary of the crystal where boundary condi- 
tions are satisfied, as long as this boundary is not 
too close to the nuclei. In fact the assumption is 
much more justifiable in the present case than in 
the case of x-ray diffraction, for in the latter 
problem the series analogous to (48) has very 
different convergence properties. In the dy- 
namical theory of x-ray diffraction, it is also 
common to assume that the surface boundary is 
a perfect plane. As in the treatment of refraction 
in the preceding section, this assumption will 
evidently be accurate only if large areas of the 
crystal surface are composed of the same plane 
of atoms in the usual crystallographic sense. 

The analogy between x-ray diffraction and 
neutron diffraction, in the ideal case we are 
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considering, extends even more deeply than might 
be supposed at first sight, as we shall see more 
clearly in the next section. 

We shall examine the equations for neutron 
diffraction further after discussing another 
method by which the problem may be treated. 


5. APPLICATION OF THE METHOD OF THE 
PSEUDO-POTENTIAL 


The results derived in the preceding sections 
can be obtained in a very straightforward 
manner with the use of the pseudo-potential, 
first employed by Fermi.? In this method the 
actual interaction potential between neutron and 
nucleus is replaced with a simple function that 
will lead to the same phase shifts as the actual 
potential for a range of energy of interest. This 
range will be taken to be that near thermal 
energies in the present problem. 

To begin with, we shall consider the case in 
which the pseudo-potential is taken to have the 
form of a simple square well of depth uh?/2m and 
breadth a. In the general case in which the 
nucleus absorbs neutrons, u will be complex. The 
relation between u and the phase shift 7 outside 
the potential well is found by requiring that the 
wave functien and its derivative be continuous 
at the boundary of the square well. The corre- 
sponding condition is'* 


sin(ua)* _—sin(ka+7) 


(49) 
ucos(ua)* k cos(ka+7) 


If (ua)*, ka, and 7 are all small compared with 
unity, both sides of this equation may be ex- 
panded. The resulting equation is 


‘n=ka*u/3 (50) 


or 
u=3n/ka*. (51) 


We note from (51) that the product ua’, which 
is proportional to the volume integral of the 
pseudo-potential is a constant in the present 
approximation if the cross sections are constant. 
Hence as long as the pseudo-potential is multi- 
plied by functions which vary very slowly over 
the potential well, as when the neutron wave- 
length is long compared with a, the square well 
may be replaced by a delta-function. In this 


approximation the pseudo-potential V; asso. 
ciated with the jth nucleus can be written in the 
form 


h? (= 2}! 
(52) 


. It should be noted again that the conditions 


under which (52) is valid are that », a, and the 
wave-length, X, satisfy the inequality relation 


n/kKa<0x. (53) 


These conditions stand in contrast with those 
for the methods employed in the previous sec- 
tions. In these cases it was only necessary that 
d be large compared with 7/k. 

It is easy to show that the relation (52) is 
valid even if the original well is not taken to be 
square, provided the phase shift is sufficiently 
small. To prove this we note that the left hand 
side of Eq. (49) is the value of R/R’ at the outer 
edge of the potential barrier, regardless of the 
form it is assumed to have, where R is r times 
the radial part of the s wave function, and R’ 
is its derivative. Now R/R’ satisfies the equation 


d/R RR" R 
=1-— =1+u(r)—, (54) 
dr\R’ R” 


in which u(r) is —2m/h? times the pseudo- 
potential. Since R/R’=r for small phase shifts, 
it follows that at the outer boundary of the 
potential (r =a) 


(=) f f u(r)dr. 


Equating this to the right-hand side of (49) we 
obtain 


1 
f side. (55) 


Thus the pseudo-potential can be taken to have 
the form (52) wherever the conditions (53) are 
satisfied. 

To summarize, we may take the Schroedinger 
equation inside the crystal to have the form 


wh 


_ 

a= 

ger 
nut 
in 

refi 
Thi 
use 
aF 
It i 

ove 

in 
kno 
stiti 
tern 

Aga 
Ay +ky = — ay, (56) 


REFRACTION OF NEUTRONS 303 


where 

= 4 


Xé(r—r,), (57) 


in which 7, is summed over all atoms in the 
lattice. It is clear that Eq. (35) is much more 
general since the only condition on its applica- 
bility is that \ be large compared with 7/k for all 
nuclei. Actually, as we shall see below, the 
results obtained with the use of Eq. (56) and the 
type of perturbation theory that was employed 
in connection with Eq. (17) for the index of 
refraction and the Laue-Bragg diffraction are the 
same as those found in the previous sections. 
This indicates that the results obtained with the 
use of the pseudo-potential are more general than 
the equations from which they are derived. 

: We shall now resolve a appearing in (57) into 
a Fourier series of the form 


a=)’ x ax exp(iK-r), (58) 


in which the ag are constants and the vectors K 
extend over the reciprocal lattice of the crystal. 
It is readily found from (57) that 


ax=4nno &(j) exp(—iK-rj), 


where mp is the density of unit cells, r; is summed 
over the atoms in the unity cell, and £(j) is 
defined by Eq. (36). 

As previously, we may now solve Eq. (56) 
with a function of the form 


ckexp(iK-r)), (60) 


in which the cx are constants, and k’ is an un- 
known vector. When (59) and (60) are sub- 
stituted in Eq. (56) and the coefficients of each 
term e®-t are equated, we obtain the following 
family of equations. 


Dx =[(k’+K)*—k*]cx. (61) 
Again there are two interesting cases to treat. 


(a) No Bragg Reflection 


In this case ¢ is larger than all other values of 
cx and Eqs. (61) may be written 


= [k’?—k? 


(62) 
axco=[(k’+K)*—k* 


The first equation leads to (41) when the relation 
(59) is employed, whereas the second reduces to 
Eq. (37) when Eq. (40) is used to express a_,(j) 
in (37) in terms of £(j). 


(b) Near Bragg Reflection 


When (k’+K’)? is near to k” for a particular 
value K’ of K, the constants cp and cg, become 
comparable. There are then two equations 
analagous to the first of (62), namely 


=[k’*—k* 
=[(k’+K’)?—k* 


These equations are identical in form with (47) 
so that the allowed values of k’ are identical in 
the two cases. Moreover, it is readily found that 
when the value (46) for a_:(j) is substituted into 
Eq. (37) for cx, the result is identical with the 
equation for cx obtained from (61), namely 


(64) 


Hence the method based on the pseudo- 
potential yields the results we obtained previ- 
ously. 

A close comparison between Eqs. (63) and the 
corresponding equations for x-ray diffraction” 
shows that the function a/4rk* can be regarded 
as the equivalent of the complex electrical polar- 
izability of the medium in which the x-rays are 
diffracted. In a region of wave-length not too close 
to an absorption band, the polarizability u(r) for 
a monatomic substance is given by the equation 


u(r) = (€*Z/mc*k*) p.(), (65) 


in which p,(r) is the electron charge distribution, 
which is normalized to unity for a single atom; 
Z is the atomic number; m is the electronic 
mass; and ¢ is the velocity of light. This is to be 
contrasted with a/4rk? which in a non-absorbing 
monatomic substance is 


(63) 


=.) (x) (66) 


where p(r) is the density of scattering centers. 
Since Ze*/mc* is of the same order of magnitude 
as (c,/4r) in regions not too close to resonance 
peaks, it follows that we may expect nearly 
quantitative agreement for the two cases. 


o 


at 
| 


This close analogy is not as valid in the 
vicinity of absorbing regions as it is the case of 
pure scattering. Absorption cross sections of the 
order of 10-* cm? or larger are very common in 
the case of x-rays, particularly when the material 
contains heavy elements. Similar large absorption 
cross sections occur for thermal neutrons only in 
crystals containing cadmium, gadolinium, or one 
of the other elements possessing a resonance peak 
in the thermal region. 


6. APPLICATION OF EQUATIONS FOR 
LAUE-BRAGG SCATTERING 


We shall now apply the results of diffraction 
theory to a few typical problems of interest in 
connection with neutron diffraction, restricting 
our attention to cases in which absorption is 
negligible, unless the contrary is stated. 


(1) Depth of Penetration 


In the conventional treatments of the dif- 
fraction of x-rays, it is customary to distinguish 
two cases, namely that of a “thin” crystal and 
that of a “thick” crystal. Only a small portion of 
the incident wave is diffracted before passing 
through the crystal in the first case, whereas the 
wave is diffracted a number of times in the 
second case. The conditions on the thickness fo 
determining whether or not the crystal is to be 
regarded as thin or thick are, respectively, that 
the quantity A defined by the equation 


A =ax|to|/2k(|yo| |yx|)# (67) 


be much smaller or much larger than unity (cf. 
Zachariasen, Eq. (3.140)). Here K is the re- 
ciprocal vector for the diffracting plane; yo is the 
cosine of the angle between the direction of the 
incident wave and the normal to the surface 
through which the wave enters the crystal; yx 
is the cosine of the angle the diffracted ray makes 
' with the same normal; and the other quantities 
are as defined previously. The two cosines are 
equal in magnitude in the case in which the 
incident wave undergoes Bragg-reflection from 
crystallographic planes that are parallel to the 
surface plane of the specimen. 

Thus the reciprocal of the coefficient of to in 
Eq. (67) provides a measure of the depth to 
which a wave penetrates before being attenuated 
to the maximum degree. This distance, which we 
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shall designate as dx is 
dx = (2k/|ax|)(| vol | (68) 
In the case of a monatomic substance this js 
dx = | 


This is closely comparable with the correspond- 
ing distance for x-rays for®® all ordinary values of 
g., as we might expect from what has been said 
previously. For example, if we take 


k=3.5-108cm™, 
o./4x = 10-*4 cm?, 
we find Vel 


d=3.7X10-* cm. 


As a more specific example let us consider the 
case of thermal neutrons diffracted by the (100) 
faces of magnesium oxide, which possesses 
sodium-chloride lattice. We shall take the neu- 
tron energy to be 3ko7T/2 (ko is Boltzmann’s 
constant, and T will be chosen as 300°K) for 
which the neutron wave number is 4.57 
cm, The most reliable values for the cross sec- 
tions appear to be o,(Mg) =4.0X10-** cm? and 
o,(0) =4.0X10-* cm*. If we assume that both 
nuclei scatter with the same phase @199 = 7.6 X 10" 
cm~*. Since (|o| | yx|)#=0.653, we obtain 


d100=7.84X 10-4 cm. 


It is to be noted that in this case doo would 
vanish to a first approximation if the phases of 
the scattering nuclei were opposite since the 
cross sections are almost identical. In fact the 
cancellation would occur for all orders of reflec- 
tion from the (100) planes since 


X(1+cosNr). (69) 


If the (100) reflections vanish, it is still possible 
to obtain reflections from the (111) planes since 


X (cos3Nx+3 cose) | (70) 
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This close correspondence between the pene- 
tration distance for the two types of radiation 
shows us that the upper limit on the size of mosaic 
blocks for which primary extinction is negligible 
for a strong x-ray reflection is about the same 
for neutrons. 


(2) The Intensity of Reflection from a Crystal 


Let us now consider the intensity of a reflected 
neutron beam under the somewhat idealized con- 
ditions we are discussing in the present treat- 
ment. To begin with, let us assume that a poly- 
chromatic beam impinges upon a perfect crystal 
and makes an angle with the plane of the crystal. 
We shall assume that this plane coincides with 
the planes responsible for diffraction so that the 
situation is that typical of Bragg reflection. In 
addition we shall assume that the number of 
neutrons per second in the beam having wave- 
length lying in the range from \ to A+) is given 
by »dd. The incident beam will possess a small 
angular divergence so that the wave-length 
which is diffracted with maximum effectiveness 
will differ for each ray of which the beam is 
composed. We shall assume that the angular 
divergence is sufficiently small that the distribu- 
tion function » does not vary appreciably over 
the range of wave-lengths which are diffracted 
most efficiently. 


According to Zachariasen (cf. Eqs. (3.152) and 


(3.156)) the integrated intensity Neg which is 
reflected from the crystallographic plane asso- 
ciated with K for a thick crystal is 


(71) 
in which \g satisfies the equation 
sinds/|K}. (72) 


Nr represents the net neutron flux per unit time 
across the entire area of the reflected beam. If 
the distribution spectrum of the incident beam 
is expressed in terms of wave number instead of 
wave-length, Eq. (71) is replaced by 


Ne ve(k)/2k sin?@z, (73) 
ve(k)dk = va (A) dX. (74) 


where 


If the incident neutron beam has a net inten- 
sity N; neutrons per second and is derived from 


a thermal source by collimation, the distribution 
function », is given by 
2N1/ 


k 2mkoT 


) exp(—h?k?/2mkoT), (75) 


in which ko is Boltzmann's constant, T is the 
neutron temperature; and the other quantities 
are as designated in previous sections. Using 
Eq. (73) we find 


Ne 
Nr k? sin?6, 2mkoT 


|ax| ye 
\RkoT 


in which ¢ is the neutron energy. If we set 
e=%koT, we find that (76), which measures the 
effective reflectivity of the crystal, is of the order 
of 10-* for a typical monatomic scatterer. Ac- 
tually the reflectivity will be this small only for 
a very perfect crystal. We may expect values 
between 10 and 100 times larger for a moderately 
imperfect crystal containing mosaic blocks which 
are tilted relative to one another. 

It is interesting to compute Ne/N; for MgO. 
We shall consider second-order diffraction from 
the 100 planes and first-order diffraction from 
111 planes. Here k=4.57X10® |ayoo| 
= |@in| =7.61X10" cm, @8100=40°45'; 
= 16°26’; and we find X10, 
(Ne/N1)in=7.10 


) exp(—A*k?/2mkoT) 
(76) 


(3) Influence of Mosaic Structure 


In a single crystal x-ray spectrometer, the 
mosaic structure seriously interferes with the 
resolving power in cases in which the highest 
detail is required, for example when it is desired 
to resolve the natural breadth of emission lines. 
The corresponding problem is probably not as 
grave for a neutron spectrometer, at least from 
present standards. For example, if we require 
an energy resolution AE in the energy range near 
E the allowed angular spread of the beam must 
be less than 


A@=(AE/E)(1/2 coté), (77) 


in which @ is the angle at which observations are 
made. When @ is about 5X10, which is a 
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typical small value that would be of interest in 
the diffraction of thermal neutrons from a crystal 
such as calcite, we find that an energy resolution 
of 5X10-*, which would be of high quality com- 
pared to present standards of resolution, can be 
achieved with an angular divergence of about 
10 radian or about 0.4 second of arc. Thus 
both ‘the degree of divergence of the neutron 
beam and the degree of variation of the orien- 
tation of mosaic blocks could be considerably 
larger for the case of neutrons than for x-rays 
for ‘‘good”’ results. 

The quantitative importance of mosaic struc- 
ture in determining the reflectivity of a crystal 
can be estimated approximately in the following 
fashion. If the crystal were perfect and were thick 
compared to dx given by Eq. (69), its reflec- 
tivity would be unity for monochromatic waves 
incident at the Bragg angle for a range of angular 
deviation w given by (cf. Zachariasen (3.151) 
and (3.155)) 


w=1.15|ax| sin26p. (78) 


This width is of the order of a second of arc for 
a typical case. In the case of magnesium oxide we 
find that when k = 4.57 X 108 cm=, w=0.875” for 
second-order reflection from the (100) planes, 
and w= 1.59” for first-order reflection from (111) 
planes. The net reflectivity for a polychromatic 
beam is given by Eqs. (71) and (73). 

The angular distribution of mosaic blocks 
relative to perfect orientation will, in the simplest 
case, be given by a distribution function of the 
form 


(79) 


in which A, is the angular deviation from the 
mean, and 7 measures the breadth of the dis- 
tribution. 7 is very small, of the order of a few 
seconds of arc or less, in a few natural crystals. 
However, in most crystals 7 is much wider, of 
the order of a minute of arc or more. 

In the case of a thick crystal (i.e., a calcite 
crystal 1 cm thick) we may expect to obtain unit 
reflectivity over a range of angle equal to 7 
instead of over the range w given by Eq. (78) 
if we assume that absorption is not large for the 
thickness of the specimen. Thus in the cases of 
interest to us, the reflectivity is equal to (76) 
times 7/w. The reflectivity then may be of the 


order of 10~‘ or larger instead of the value given 
by Eq. (76). 

In the corresponding problem for X-rays, the 
absorption of the incident beam is usually 


great that the simple reasoning of the previous 


paragraph cannot be employed. It is readily 
shown that when the absorption is great the 
quantity multiplying (76) is (cf. Zachariasen 
(4.40)) 
(o./ koa, (80) 
provided (Zachariasen (4.34)) 
2x? m 1 af 


(81) 


in which 7 is expressed in radians. The condition 
(81) is generally not met for neutrons in the 
cases of interest. For example, 7 is of the order 
of 3X10-*, whereas the expression on the left- 
hand side is of the order of 10-*. The approxima- 
tion (81) becomes valid if o is about 10* times 
o,. However (80) is then near unity. 


7. ADDITIONAL INFLUENCE OF DISORDERING 
UPON THE INTENSITY OF BRAGG 
REFLECTION 


We have seen in the previous section that the 
angular disorientation of the mosaic blocks of 
which real single crystals are composed increases 
the reflectivity of a specimen for polychromatic 


‘ radiation by increasing the range of angles of 


incidence for which a monochromatic com- 
ponent can be reflected. Most other types of 
crystalline disorder decrease the reflectivity. 
Two disorders of this type which should be 
mentioned are the disordering arising from the 
thermal (and zero point) oscillations of the atoms 
and the disordering which can occur when the 
atoms are polyisotopic. 

The influence of atomic oscillation has been 
discussed in detail by Weinstock” for a mon- 
atomic lattice in the case in which primary ex- 
tinction is negligible, that is, in the case in which 
each mosaic is so small that only a small fraction 
of the monochromatic radiation which satisfies 
the Bragg condition is scattered. In this case the 
reflectivity is reduced from the value for a 
crystal in which the atoms are at rest by a factor 


(82) 
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_where )o(=)] (83) 
with 
Q(z) = J an 


Here K is the vector in the reciprocal lattice 
associated with the reflecting plane; @ is the 
Debye characteristic temperature; ko is Boltz- 
mann’s constant; and M is the mass of the atoms 
of which the lattice is composed. It is readily 
seen that W is of the order of m/M for thermal 
neutrons when 7 is near @ and @ is near room 
temperature, as is the case for the crystals with 
comparatively large elastic constants that would 
normally be used in a neutron spectrometer. m is 
the neutron mass. Thus the factor (82) is usually 
very close to unity and will not influence the 
conclusions drawn previously in an appreciable 
way. - 

We shall discuss the influence of isotopes for a 
monatomic lattice since the results can readily 
be generalized. Let us consider a small crystal 
upon which is incident a wave possessing the 
wave number vector k. The quantity which de- 
termines the intensity of reflection in the direc- 
tion for which the wave number vector is k’ is 


Lif; (84) 


in which r; is the position of the jth atom in the 
lattice, and f; is proportional to +(¢,*)+. If the 
element is monoisotopic, the factor f; is inde- 
pendent of 7 and may be removed from the 
summation. The square of the absolute value of 
the resulting summation may then be approxi- 
mated by a sum of delta-functions with constant 
multipliers which are centered about the points 
for which k—k’=K, where as before the K are 
the reciprocal vectors. On the other hand if two 
or more isotopes are present, and if the values of 
fi for these isotopes are different, the value of 
(84) in a given direction for which k—k’ =K is 


Ls (85) 


where N is the number of atoms, and f is the 
average value of f;. This average is taken in such 
a way that each isotope is weighted with its 
abundance and due recognition is taken of the 
sign before (c,’)*. Thus even when the values of 


o,/ are nearly identical for the isotopes the re- 
flectivity may be decreased, provided the major 
isotopes are present in comparable abundance and 
have opposite phases. Conversely the isotopic 
mixture will be unimportant if all isotopes have 
nearly equal cross sections and the same phases. 

It is readily shown that the distribution of 
the scattered radiation is described by a set of 
delta-functions peaked at positions for which 
k—k’=K upon which is superimposed a diffuse 
isotopic component of scattered radiation. The 
diffuse component is isotopic only when the 
isotopes are distributed at random among the 
various lattice sites. To demonstrate this, we set 


fi=I+ (86) 


where >>; 6f;=0. We may then write (84) in the 
form 


Xs expli(k—k’)- 

+2; 6f;exp[i(k—k’)-r;]. (87) 
The first sum is sharply peaked at k—k’=K, 
as when the 4f; are zero. Hence any diffuse 
scattering must be described by the second 
term. The square of the absolute value of this 
term is 

5fj6f,* exp[i(k—k’) -(r;—ry) ] 

which may also be written in the form 


(88) 


where r,=1j—1,,,. The summation for fixed s 
may be separated from that over j, that is the 
summation may be written as 


Le exp[i(k —k’) 5f (89) 


When the isotopes are randomly distributed, the 
sum over j vanishes if s is not zero since for each 
particular value of éf; that occurs f;,, will range 
over all values at random and the resultant is 
zero. When s vanishes the second sum is a finite 
constant which is, of course, independent of s 
and may be removed entirely to the left of the 
summation over s. In this case r,=0 so that the 
expression for the scattered radiation is inde- 
pendent of (k—k’). Thus in addition to the 
Laue-Bragg scattering there is an isotropically 
distributed diffuse component. 
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If the isotopes were correlated, >>; dfi5fi,.* 
would be a function of s which is different from 


zero when s#0 and diffuse, but anisotropic scat- 


tering would occur. This scattering would re- 
semble closely the type of diffuse x-ray scattering 
observed for liquids. There is, of course, no 
reason for expecting a correlation between the 
positions of isotopes in the heavier solids. 


SUMMARY 


1. The equations for the elastic scattering of 
neutrons by a single crystal are examined in 
detail in the approximation in which inelastic 
scattering is neglected. This procedure is equiva- 
lent to treating the case of a perfectly ordered 
crystal of rigidly fixed atoms by methods 
analogous to those employed in the dynamical 
theory of x-ray scattering. Three methods of 
deriving the results are employed: The first, 
which is the most rigorous method takes account 


of the singularity produced on the wave function 


in the region near but just outside the nucleus 
in the general case in which the nucleus both 
scatters and absorbs. This method is accurate as 
long as the square of the wave-length is large 
compared with the cross sections. The second 
method is based on the use of the cellular ap- 
proximation that has found wide use in the 
treatment of electron wave functions in solids. 
It is applied only when the absorption cross 
section vanishes, but is then as accurate as the 
first method. The third method is based on the 
use of a pseudo-potential, which replaces the 
actual potential between neutron and nucleus 
in the Schroedinger equation. This method, 
which was first used by Fermi, can be justified 
only when the wave-length is large compared 
with the diameter of the pseudo-potential well, 
which must in turn be large compared with the 
“displacement distance.” 

2. It is shown that the index of refraction n 
and the extinction coefficient mx are given by the 
following equations when the wave-length and 
direction of incidence of the neutrons are not 
such that the Bragg conditions are satisfied: 


fos) 
n=1+m j;+—{ ——F? ) 
02 4r 
nk=MNy 


These are generalizations of a similar set of 
equations derived by Fermi. Here mp is the 
number of unit cells per unit volume of the 
crystal; k is the wave number, which is equal to 
where is the wave-length; is the scat- 
tering cross section of the jth atom in the unit 
cell; and go,’ is the corresponding absorption 
cross section. The signs which appear before the 
radicals in the first equation are related to the 
sign of the shift in phase of the wave scattered 
by the corresponding atom. These equations 
may be used for highly disordered materials jn 
the case in which the wave-length is so large 
that the material is homogeneous over a distance 
comparable to the wave-length. In this case My 
is to be taken as the number of chemical mole. 
cules per unit volume, and the summations 
extend over the constituents of the chemical 
molecule. It is interesting to note that the second 
equation can be written in the form 


in which L, is the mean free path for absorption 
in the material. 

3. A crystal will reflect neutrons which strike 
it at grazing incidence provided the index of 
refraction is less than unity. This is to be ex- 
pected from the general laws of optics and has 
been pointed out previously by Fermi in con- 
nection with is derivation of the equations dis- 
cussed in paragraph 2. The reflectivity for a 
monochromatic beam of neutron is given by the 
equation 


—(1 
1+(1+2y/¢?)! 


in which ¢ is the angle that the beam makes 
with the reflecting surface and ~4=1—n. Ab- 
sorption influences the reflectivity only when 
there is a resonance level in the thermal region. 
The depth of penetration of the ray reflected 
at grazing incidence is given by the equation 


This distance is of the order of magnitude 100 
for a typical case. 

The experiments carried out by Fermi, Zinn, 
Marshall, Borst, and Wollan seem to show that 
the index of refraction is usually less than unity. 
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It is pointed out that this observation may be 
no more surprising than the fact that most 
transparent optical media have an index greater 
than unity in the visible portion of the spectrum. 

4. It is found that there is a very close cor- 
relation between the diffraction of neutrons and 
of x-rays which are polarized in the direction 
pormal to the plane of incidence. This correlation 
is almost quantitative in the non-absorbing cases 
because the cross section for Thomson scattering, 
which is related to the classical electron radius, 
is almost equal to the typical cross section for 
nuclear scattering. The absorption cross section 
for neutrons is comparable to the corresponding 


typical cross section for x-rays only in the special. 


case in which one of the nuclei present possesses 
a resonance level in the thermal region of the 
energy spectrum, that is, when elements like 
cadmium or gadolinium are present. Only in 
cases of this type does the absorption of neutrons 
affect the properties of a material for diffracting 
neutrons. 

-5. The depth of penetration of a diffracted 
wave which satisfies the Laue-Bragg condition is 
determined for a perfect crystal. The result is 


d=2k/|ax| 


in which 6g is the Bragg angle and 


Xexp(—iK-r,). 


Here all quantities are as in paragraph 2 above, 
and K is the reciprocal vector of the lattice that 
is to be associated with the crystallographic 
plane from which the incident beam is reflected 
in the Bragg sense. 

a plays much the same role in neutron dif- 
fraction as the x-ray structure factor does in 
x-ray diffraction. 

In the case of a non-absorbing monatomic sub- 
stance this equation simplifies to 


sin@g 
ANo 


This distance is of the order of 10~ in typical 
cases, just as in the case of strongly reflected 
x-ray beams. 


6. The reflectivity R of a perfect single crystal 
for a stream of thermal neutrons is given by the 
equation 


2 
~) 
k? \koT 


Here ko is Boltzmann’s constant, and ¢ is the 
energy of the neutrons which undergo Bragg 
reflection for the geometrical arrangement em- 
ployed. This reflectivity is of the order of 0.001 
percent when the angle of incidence is adjusted 
so that the neutrons at the peak of the thermal 
distribution function in the incident beam are 
diffracted. 

7. It is pointed out that the reflectivity will 
be much larger than the equations given in 
paragraph 6 indicate if the crystal possesses well- 
developed disorientation of the type in which 


the mosaic blocks are tilted relative to one. 


another. If the mean spread in angular disorien- 
tation is 7, the expression for the reflectivity 
given in the preceding paragraph should be mul- 
tiplied by a factor 7/w where w is the rocking 
angle for which high reflectivity is obtained 
when a highly collimated monochromatic beam is 
incident on a perfect crystal near the Bragg angle. 
w is of the order of several seconds of arc. This 
multiplicative factor is to be used only when 
n/w>>1. Thus we see that by selecting a crystal 
with the proper mosaic structure, it is possible 
to enhance the reflectivity by a factor of one or 
two decades. 

8. If it is desired to maintain a resolving 
power that is considered ‘“‘good’’ from present 
standards, that is, an energy resolution of the 
order of 5X10 ev for 1 ev neutrons, the mean 
angular disorientation caused by the mosaic 
structure must be no more than a few minutes 
of arc. The angular aperture of the entrance and 
exit slits of the system should also be no larger 
than this if the same resolving power is specified. 

9. The decrease in intensity of Bragg reflection 
arising from the thermal oscillations is discussed 
on the basis of Weinstock’s calculations. It is 
shown that this decrease is not important for 
crystals which possess high elastic constants, as 


those used in a neutron spectrometer probably * 


would. 
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10. It is shown that the reflectivity may be 
decreased if the atoms in the crystal are polyiso- 
topic and if the isotopes for a given element are 
present in comparable abundance, possess nearly 
equal scattering cross sections and opposite 
phases. More explicitly the effective scattering 
strength of a given element is proportional to 


i, 
in which g; is the fractional abundance of the ith 
isotope, and a,‘ is its scattering cross section. 


The contribution from isotopes with opposite 
phases tend to cancel one another. The presence 


HARRIS 


of isotopes may also add a diffuse background to 
the inelastically scattered radiation. 

11. It should be emphasized that the present 
paper does not contain a discussion of the go. 
herent inelastic scattering of neutrons from 
crystals. This type of scattering actually can be 
very important in practical cases and can alter 
some of the qualitative viewpoints which were 
obtained. 

We are indebted to our colleagues, L, Borst, 
S. M. Dancoff, E. P. Wigner, E. O. Wollan, W. H. 
Zachariasen, and W. H. Zinn for many interest. 
ing discussions of this topic. 


PHYSICAL REVIEW VOLUME 71, NUMBER 5 MARCH 1, 1947 


A Proposed Focusing Cosmic-Ray Telescope 


W. T. Harris 
United States Time Corporation, Waterbury, Connecticut 
(Received October 31, 1946) 


A magnetic lens which will focus charged particles entering its aperture parallel to the axis 
can be constructed in the form of a toroidal winding. For an air or iron core toroid, the cross 
section of the winding is parabolic. If partial iron filling is used, trapezoidal or rectangular 
cross sections may be employed to produce sharp focusing. The focal length is proportional to 
the particle momentum, and hence a telescope of this type is also a spectrograph. By using 
coincidence counters and anticoincidence guard counters for observations, a collimating tube 
effect can be obtained, permitting full use of the properties of the telescope lens without a 
background due to particles which do not pass through the lens. There are some advantages 
to an installation in a deep shaft in the earth having two lenses. With such a device charged and 
uncharged components can be studied separately. 


INTRODUCTION 
A° a consequence of Maxwell’s equation 


curl H=4zi, 


it is clear that the magnetic field strength inside 
a long straight solenoid is uniform over the cross 
section of the solenoid, regardless of the shape of 
the cross section. If the solenoid is formed into 
a toroid, H varies inside the cross section as 1/r, 
regardless of the shape of the cross section of the 
toroidal winding. 
If a toroidal coil is to serve as a lens to focus 
»? parallel beam of particles of like energies and 
e/m through the point F (Fig. 1), it is necessary 
that the radial momenta imparted to the par- 


ticles increase directly as r. This can be accom- 
plished in either of two ways. 

(1) Air core toroid. Since in an air core toroid 
the field varies as 1/r, the width of the cross 
section must increase as r?, in order that the time 
in which a particle is accelerated radially may 
compensate the decrease in field strength and 
yield a radial impulse which is proportional to r. 
This leads to a toroidal lens having the parabolic 
shape shown in Fig. 1. 

(2) Toroidal lens containing iron. If iron is 
distributed inside the toroid, and the iron is 
regarded as being transparent to the particles, 
several choices of lens design are possible. In this 
case we must deal with the continuous magnetic 
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flux B in the magnetic circuit, rather than the 
magnetic field H.' 

(a) The interior of the toroidal coil may be 
uniformly filled with iron. The resulting flux dis- 
tribution is the same as for the air core case (1) 
above. (The situation is actually the same as if 
the permeability of all space were changed by a 
uniform amount.) Thus the only advantage of 
the use of iron in this way is the large saving in 
power required to operate the lens, due to the 
smaller current necessary to obtain the desired 
flux. 

This type of lens would not be very useful 
when small magnetic fields are required. In addi- 
tion to the collision effects in the iron, poor 
reproducibility of results might constitute a 
further disadvantage. Since the magnetic circuit 
would be completely closed, the demagnetizing 
force would be low, and the slight coercive force 


MAGNETIC FIELD INTO PAPER 
CONFINED TO TORQIOAL SPACE 
(NSIDE WINDING AS SHOWN- 


Fic. 1. For an air core lens or a homogeneously iron filled 
lens, the toroidal magnetic field has a parabolic cross 
section. 


? Recently some doubt has been raised regarding the 
equivalence of the interaction of charged particles with a 
given magnetic flux inside iron and in space. If these 

posals had any real basis the two partially iron filled 
orms of telescope lenses described below mene « f be unsatis- 


factory. See G. Wannier, Phys. Rev. 67, 364 (1945) and 


D. L. Webster, Phys. Rev. 70, 446 (1946). 


of the material (ca. 0.1 oersted) would be suf- 
ficient to prevent the flux from being a unique 
function of the magnetizing field at low flux 
densities. However, complete iron filling will be 
necessary for particles of very high momentum, 
such as 10'°-volt protons. 

(b) The design is more flexible, and the above 
objections are overcome, if the toroid is only 
partially filled with iron. If, for example, the 
volume inside the winding is half filled with iron 
for which u>1 (ca. 3000), the interior can be 
regarded as having a permeability of 2 for the 
purpose of calculating total flux. We are now 
privileged to choose to make the flux density 
uniform inside the toroid, which requires that 
the toroid have a trapezoidal cross section for 
focusing. The iron and air volumes would alter- 
nate, the iron pieces having the shape of trape- 
zoidal wedges, and the air gaps would be 
trapezoids of uniform width, as shown in Fig. 2. 
A large number of relatively thin wedges should 
be used to reduce fringing flux effects. 

(c) A toroid of rectangular section can be used 
if the iron is so disposed that the flux density 
increases directly as r. This result can be achieved 
by increasing the slope of the wedges, so that the 
air gaps become narrower in proportion as r is 
increased. See Fig. 3. 


Fic. 2. Toroidal magnetic lens partially iron filled so as 
to produce a constant flux{density. In this case the lens 
thickness must increase directly as r. 
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Fic. 3. Partially iron filled lens in which the magnetic flux 
increases as 


In order to test the practicability of the 
toroidal lens the magnetic field strength require- 
ments are of primary interest. As an example, we 
assume an f:5 lens for 10%-volt protons 50 cm 
thick at the outer periphery and half filled with 
iron. (See Fig. 4.) Taking 8 as 0.44, the time to 
traverse the magnetic field at the outer edge of 
the lens is given by 50=0.44ct, or t=3.8X10-° 
sec. The radial momentum, p,, imparted to the 
particle is (Bev/c)t. If we set this equal to 0.1 of 
the momentum of the particle Hpe/c, we have 
3.8X10-°Bv=0.1Hp. For 108-volt protons Hp is 
approximately 1.4x10°. Hence the magnetic 
flux required is 


B=1.4X105/(3.8X10-X 1.3 10") = 2800 gauss. 


Hence the requirements are reasonable, although 
at energies between 10° and 10!° ev it may become 
difficult to make the focal length for heavy par- 
ticles short enough to preserve the usefulness of 
the device. 


CONSTRUCTION 


As visualized at present, the magnetic lens 
should consist of a large wound doughnut, pre- 
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OF COUNTER 


Fic. 4. Sketch illustrating computed example. Counter 
arrangement to prevent detection of particles which do 
not pass pews + the lens is shown. Instead of cylindrical 
counters, parallel disks with a few fine wires in the median 
plane between them would be more suitable. 


ferably with the windings open and arranged 
for forced air cooling. Square section copper bars 
would be employed for the winding, and the size 
of the center hole would be determined by the 
construction. The windings might be single layer, 
close wound at the inner limit, and spaced at the 
outer periphery, or they might be multilayer at 
the inner limit with a smaller number of layers 
at the outer limit. The use of iron in the lens has 
previously been discussed. For studying low 
energy particles an air core lens would be desir- 
able, but for studying very high energy particles 
iron would be desirable. 

For most purposes, lenses probably should be 
rather large. A }- to 1-meter radius for the inner 
hole, and a 3-meter radius to the outer periphery 
is suggested. This would have a useful focal 


range of perhaps 5 to 50 meters. The projected : 


area of the active lens would be of the order of 
10 square meters. This would give an expected 
concentration of a homogeneous component 
coming from above of approximately 5000, if 
good focusing and a coincidence counter bank of 
20 cm? active area are assumed. (It should be 
noted, however, that the earth’s magnetic field 
will displace and blur the focal spot. A rather 
clumsy Helmholtz type coil arrangement can be 
used to approximately neutralize the earth’s 
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field in the telescope “tube”’ and remove this dif- 
ficulty.) This concentration, the high directivity, 
and the chromatic abberation which renders the 
device a spectrometer are the particular advan- 
tages of the device. 

The ‘power required to produce a uniform 
magnetic flux of 2800 gauss in a toroid of 3 meters 
outer radius and center hole of 1 meter radius, 
the toroid being half filled with iron, and single 
layer wound with 1” square copper bars (or an 
equivalent volume of copper) as suggested above 
is 970 kw. If the toroid were completely filled 
with iron the power required would be less than 
1 kw. Increasing the volume of copper decreases 
the power consumption in direct proportion. 

It is well to point out that the above discussion 
has been confined to particles which approach 
the lens parallel to its axis. In order to study 
particles which approach at only a small angle 
to the axis in an experimental set-up it is neces- 
sary to eliminate detection of particles which 
come from other directions. This can be accom- 
plished by two different methods. (1) Banks of 
coincidence counters which count only particles 
passing through the aperture of the lens can be 
employed, together with off-axis anticoincidence 
counters to eliminate counts due to showers. For 
the coincidence counters parallel plate pancake 
counters with an open grid of wires in the central 
plane are suggested. In this way a collimating 
tube effect is obtained, and full use of the direc- 
tivity, dispersion, and focusing effect of the lens 
may be utilized. Experiments using absorbers 
between the counters would assist in determining 
the masses of the particles. 

. (2) A second collimation method which might 
be feasible would utilize a deep mine shaft ‘in 
which to install the telescope. Some such sites 
may be available at moderate elevations above 
sea level. 

The mine shaft type of installation, of course, 
could not be trained. However, it apparently has 
advantages. For example, if the shaft were suf- 
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Fic. 5. Sketch of a proposed installation in a deep shaft 
employing two lenses. Lens A serves as a charged particle 
spectrograph and also deflects cha particles out of the 
“tube” so that the absorber B and lens C can be used to 


study uncharged components. 


_ ficiently deep, two telescope lenses could be 


installed, one above the other. (See Fig. 5.) The 
upper lens could be used to study charged par- 
ticles and as a filter to deflect charged particles 
out of the aperture of the lower lens, and hence 
permit more clearcut use of the lower telescope, 
with an absorber above the lens, to study un- 
charged components. 

If there were any reason to believe that an un- 
charged primary cosmic-ray component existed, 
this feature might be particularly important, 
since it would permit the study of discrete 
sources, if such sources exist, without the com- 
plicating effects of magnetic fields in space. In 
any case, however, the telescope with two lenses 
should be particularly useful in studying un- 
charged secondary radiation. 
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PHYSICAL REVIEW 


N a recent experiment Conversi, Pancini, and 


positive and negative mesotrons coming to rest 
in iron or in graphite. They find that in iron the 
disintegration electrons are observed only for 
positive mesotrons. This was indeed to be ex- 
pected? because negative mesotrons after being 
slowed down can approach the nuclei and dis- 
appear by nuclear interactions. If, on the other 
hand, graphite is used for stopping the meso- 
trons, delayed disintegration electrons are ob- 
served to be about equally abundant for positive 
and negative mesotrons. This is in sharp dis- 
agreement with current expectations and seems 
to indicate that the interaction of mesotrons 
with nucleons according to the conventional 
schemes is many orders of magnitude weaker 
than usually assumed. The disappearance of a 
negative mesotron can be analyzed into a process 
of approach of the mesotron to the nucleus and 
the process of capture by short range interaction 
of the mesotron and the nucleons. 

The slowing down of mesotrons to an energy 
of about 2000 ev takes place according to the 
conventional theory. In estimating the energy 
loss for lower energies we have considered energy 
exchange with electrons and radiation. We con- 
sider the electrons as a degenerate gas with a 
maximum velocity vo and assume that the 
velocity V of the mesotron is small compared to 
vo. Then the energy loss to the electrons per unit 
time is of the order of magnitude e*m?7T/(h®z). 
Here m and yu are the masses of the electrons and 
the mesotrons, respectively, and T is the kinetic 
energy of the mesotron. This formula allows 
losses of energy even when the total energy is 
negative (mesotron bound to an atom), and is 


1M. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 
71, 209 (1947). 
2S. Tomonaga and G. Araki, Phys. Rev. 58, 90 (1940). 
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Piccioni! observed separately the behavior of 
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The Decay of Negative Mesotrons in Matter 


E. Fermi, E. 
University of Chicago, Chicago, Illinois 
AND 


V. WEISSKOPF 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received February 7, 1947) 


valid as long as the mesotron moves outside the 
K orbit. At closer distances the formula will be 
somewhat modified and at the lowest energies 
loss by radiation will predominate. The mesotron 
reaches its lowest orbit around the nucleus in 
most solids in not more than 10-” second. This 
orbit is 200 times smaller than the radius for the 
K shell, which is for carbon about 10 times the 
nuclear radius and for iron about twice the nu- 
clear radius. After reaching this orbit the meso- 
tron can be found within the nucleus with a 
probability of 1/1000 in the case of carbon and 
a probability 75 in the case of iron. 

According to the conventional mesotron theo- 
ries, one will have to assume that the capture 
now proceeds according to one of the following 
schemes: 


P+yp-=N+hy (1) 
X+p-=N+Y. (2) 


Here P and N stand for proton and neutron, 
signifies the mass of the mesotron, hy is a light 
quantum, and X and Y stand for initial and 
residual nuclei in the capture process. The first 
calculation of these processes for a special form 
of mesotron interaction is due to Kobayasi arid 
Okayama, and Sakata and Tanikawa.* The re- 
sults depend to some extent upon the spin of 
the mesotron and the form of the interaction 
assumed. For example, in the case of pseudo- 
scalar mesotrons with an interaction energy 
given by S (p*oy) grad ¢; is 
the wave function of the nucleons, ¢; of the meso- 
trons, » is the mesotron mass, the index i refers 
to the charge, and r is the isotopic spin operator), 
one obtains for the time of capture by process (1) 
for a mesotron already captured in its lowest 


*Kobayasi and Okayama, Proc. Phys.-Math. Soc. 
th. Soc. Japan 21, 58 (1939). 


aes 21, 1 (1939); Sakata and Tanikawa, Proc. Phys.- 
a 
14 
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orbit 10-!* and 10~*° second in carbon and iron, 
respectively. Process (2) is likely to lead to 10 
times shorter lives. This is negligible compared 
to the life of a negative mesotron of 2X10-° 
second. . 

The experimental result! leads to the con-" 
clusion that the time of capture from the lowest 
orbit of carbon is not less than the time of natural 
decay, that is, about 10~® second. This is in 
disagreement with the previous estimate by a 
factor of about 10". Changes in the spin of the 
mesotron or the interaction form may reduce 
this disagreement to 

If the experimental results are correct, they 
would necessitate a very drastic change in the 


forms of mesotron interactions. The result is 
significant also for the production of single 
mesotrons by artificial sources. Indeed the crea- 
tion of a mesotron by x-rays or fast protons 
is the reverse of processes (1) and (2). If the 
interaction according to these two processes is 
much weaker than expected, one would conclude 
the same for the reverse processes. Thus one 
might be in doubt as to whether one can pro- 
duce abundant numbers of artificial mesotrons 
with bombardment-energies only a little above 
the threshold for single-mesotron production. 
Predictions concerning the creation of mesotron 
pairs by electromagnetic radiation are, of course, 
not affected by these arguments. 
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The Lateral Extension of Auger Showers 


D. V. Skospe.tzyn, G. T. ZATSEPIN, AND V. V. MILLER 
P. N. Lebedev Physical Institute of the Academy of Sciences of U.S.S.R., Moscow, U.S.S.R. 


(Received January 7, 1947) 


N 1939 P. Auger and co-workers discovered the sults for much larger distances. We applied a 
large atmospheric showers of cosmic rays different method of observation the essence of 
which cause the coincidences in two Geiger- which will be clear from Fig. 1. In this figure 
Miiller counters separated by a distance of sev- 1, 2, 3, 4 are four trays of Geiger-Miiller counters; 
eral dozens of meters. They also obtained some I and II are circuits which record the double co- 
evidence that the double coincidences C, can be incidences (1, 2) and (3, 4); III is a circuit which 
observed at much larger separations, as much as__ records coincident output pulses from I and II 
300 meters (the coincidence rate due to showers, and thus records fourfold coincidences (1, 2 
being in this case, however, of the same order +3, 4). 
as that of the random coincidences).' The effective area of each counter tray was 
As a result of numerous theoretical investiga- 1840 cm*. The apparatus was placed in light 
tions of the Auger showers it was concluded that veneer cabins. The pulses were sent through high 
the latter are of the usual cascade type but of frequency cables. 
ultra-high energy. In particular Auger’s curve Auger recorded the simultaneous passage of 
C,=f(D) seems on the whole to be in good accord two particles while we registered two simul- 
with the assumption of the cascade nature of the taneous pairs of particles. This considerably 
showers and in agreement* with the predictions decreased the number of random coincidences. 
of the cascade theory. However the point at 
D=300 m apparently upsets somewhat this 
harmony. 
Auger’s method could not yield reliable re- 


' Auger, Maze, and Robley, Comptes rendus (Paris) 
208, 1641 (1939), 
* D. V. Skobeltzyn, Comptes rendus USSR 37, 14.(1942). 


Fic. 1. The arrangement of four counter trays (1, 2, 3, 4) 
and the coincidence circuit. 
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On the other hand the large effective area of our 
counter trays permitted us to obtain much higher 
counting rate (compared to those of Auger). 
For a resolving time in circuit III of r~4.10-* 
sec. the number of random coincidences at 
D=300 meters was less than 0.7 percent of the 
number of true coincidences; at D=600 m the 
shower counting rate was 0.6 per hour while the 
random coincidences counting rate was 0.03 per 
hour. 

The results of the measurements are illustrated 
by the curve in Fig. 2. We thus have here re- 
liable evidence in favor of the existence of four- 
fold coincidences at D =600 m, and there is also 
evidence that these coincidences can be ob- 
served at D= 1000 meters. (Because of a shortage 
of time we obtained only three coincidences at 
this distance.) The. method applied by us per- 
mits one to carry out observations at even 
larger separations. 

The mean radius of Auger showers as com- 
puted in a number of papers seems to be about 
100 meters. The existence of fourfold coincidences 
at 600 meters (i.e., at separations exceeding by 
6 times the shower radius) sharply contradicts 
these calculations. 

The observed discrepancy may be estimated if 
one makes some definite assumption as to the 
dependence of the shower density p on the dis- 
tance r from the shower axis, for distances 
r>R. If, for example, one supposes that 
p~exp[—r/R]/r for r>R* (this is practically 


Fic. 2. The logarithm of the coincidence rate C, versus 
distance D in meters. The figures on the curve are the rate 
C, of fourfold coincidences per hour. 


* The assumption p~ exp[—r/R] probably involves an 
overestimation (and a very considerable one) of the density 
e at — distances r. This seems to be certainly the case 

or r 


equivalent to the assumption p~exp[ —r/R] tin 
r> R) than the ratio C4(D = 100 m)/C,(D = 600 m) 
may easily be calculated and compared with the 
experimental values. As a matter of fact on such 
assumption the function C,(R, D) for fourfold 
coincidences (at D>R) coincides with that 
which one of us found for double coincidences 
assuming p~exp[ —r/R]/r. 

The formula for C, can be obtained by sub. 
stituting in the formula previously deduced for 
C2(R, D) half the shower radius, i.e., 


C.(R, D) =C2(3R, D) 
or 


Co ff 


*texp[ —2x/R] 
= k’ 
x+D cos¢ 


 exp[—é]dt 
=2rk’ 
— (2D/R)*}* 


where k’ is a factor which depends on a certain 
integral of E. 


 p2=p0 exp(—r2/R)/r3 


are the densities of the shower particles at the 
points at which the counter trays are located, 
i.e., at distances r; and r2 from the element of 
area dS of the shower axis. F(EZ)dE is the number 
of showers of primary energy E per unit area. 
Integral (1) is the tabulated Macdonald’s func- 
tion (or 44H '(2Di/R) in which Hy)’ is the Bessel 
function of the third kind). 

The observed curve and function (1) are 
plotted in Figs. 2 and 3 on a semi-logarithmic 
scale. 

The curves in Fig. 3 were made to coincide 
at D=100 m. 

At D=600 the discrepancy reaches a value of 
about 500 times. If R=200 m this discrepancy 


** It is assumed that pio and pee<1. Such a condition 
is fulfilled for D>R if the energy spectrum is of the form 
F(E)~E-*4, This conclusion may be however invalidated 
if the number of primary particles in the region of E=10" 
— 10" ev is much greater than that which ‘can be deduced 
from this relation. 
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would be much smaller. In order to fit the part 
of the curve lying between D=300 m and D 
=600 m with the calculated curve one must 
assume R=350 m. This assumption, however, 
would lead to other difficulties and in particular 
would violate the agreement between the ex- 
perimental and theoretical curve C2(R, D) for 
D<100 m (Auger’s curve is in satisfactory 
agreement with C:(R, D) for D lying between 4 
and 100 m)?. 

This result seems to raise the question as to 
the accuracy of the usually accepted conception 
concerning the production of the Auger showers. 

If, for example, the shower-producing elec- 
trons are themselves produced in the atmosphere 
of the earth as a result of some hitherto unknown 
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process then it may be possible that at the very 
birth of the shower several particles (electrons 
or photons) of ultra-high energy are created 
simultaneously in the atmosphere, their axes 
diverging to such an extent as to be at consider- 
able distances from each other (of the order of 
several hundreds of meters) after passage through 
the atmosphere, even if the initial angle of di- 
vergence is of the order of one degree.t The co- 
incidences at large distances could then be caused 


t This angle however is much greater than mc*/E. 


Fic. 3. Logarithm of fourfold coincidence rate C, versus 
distance (experimental and theoretical curves). The figures 
on the curve are counting rates Cy per hour. 


by different but “correlated” showers. A direct 
experimental check of this point is being planned 
at present. 

This investigation was carried out in 1946 at 
the Pamir mountains at an altitude of 3860 m 
above sea level. 
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PHYSICAL REVIEW 


Letters to the Editor 


UBLICATION of brief reports of important dis- 

coveries in physics may be secured by addressing them 
to this department. The closing date for this department is, 
for the issue of the Ist of the month, the 8th of the preceding 
month and for the issue of the 15th, the 23rd of the preceding 
month. No proof will be sent to the authors. The Board of 
Editors does not hold itself responsible for the opinions ex- 
pressed by the correspondents. Communications should not 
exceed 600 words in length. 


The Origin of Excess Ultrasonic 
Absorption in Water 


LEONARD HALL 
Brown University, Providence, Rhode Island 
February 1, 1947 


HE origin of the excess ultrasonic absorption in many 

liquids very likely lies in a lag in the rearrangement 

of the molecules during a compression. This note sketches 

the theory upon which an approximate calculation of ab- 

sorption in water has been made which agrees with experi- 
mental values within the accuracy of such data. 

Under acoustic compression some of the water molecules 
break their intermolecular hydrogen bonds and jump to 
sites of closer packing. Characteristic of such transitions 
are a decrease of volume and a rate-determining activation 
energy. The activation energy should be close to that for 
shear viscous flow, while the configuration of closer pack- 
ing should be approximately that of closest possible pack- 
ing, hexagonal or cubic. 

The acoustic wave perturbs the transition rate to and 
from closer packing through the work performed by or 
against it in the volume change. Assuming the possibility 
of representation by two energy states, the unperturbed 
reaction or transition equation, —dN,/dt=kiNi—k:N2, 
will have its rate constants, k; and ke, and populations, 
N; and N2, perturbed by the acoustic radiation to k; 
etc. The excess 
pressure is p, while N,+N:2 is Avogadro’s number. From 
the solution of this perturbed reaction equation for (8Ni/ 
8p)o one obtains for the compressibility 


1 
Ni) ap 1+cos(AF/RT) (1/ks") 
Br 


where AF is the difference in free energy between the two 
states; and AV/N is the volume change per molecule- 
transition. 

This mechanism is equivalent hydrodynamically to a 
relaxational or complex isothermal compressibility. The 
rapid compressibility, 6;, arising from a uniform decrease 
in the intermolecular distances (aside from molecular re- 
arrangement) such as occurs in crystals, may be estimated 
from the Bernal-Fowler model of water or the compressi- 
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bility of ice. Since the relaxational compressibility, g, ig 
also =fo—f;, where Bo=the static value, 6; can aid in 
the evaluation of quantities occurring in the derived ex. 
pression for 
Neglecting small effects, the absorption due to re 
arrangement lag is given by 
2a/v* = -B,/cBo 
2a being the amplitude absorption coefficient (per centj- 
meter), », the frequency, and c, the velocity of sound, 
From the rate expression, k; = (kT/h) exp—[AF*/RT}, 
given by Eyring in absolute reaction rate theory, the re. 
laxation time, r= (1/k;°+1/k2°), may be determined and js 


where AF* is the free energy of activation for flow. 

Up to the present, experimental error in absorption 
measurements in water has been estimated at 30 percent, 
The values computed thus far agree most closely with the 
observed values of Fox and Rock recently published in 
The Physical Review. A complete account is in preparation, 

The present approach should be applicable to other 
liquids. 


_ A Note on Rotational Motion in 
Relativity Theory 


E. L. Hit 


Department of Physics, University of Minnesota 
Mineapolis, Minnesota 


February 3, 1947 


N a recent paper Rosen! has examined the rotating disk 
problem in relativity theory from the point of view of 
Born’s 4-dimensional definition of rigid body motion. This 
gives a welcome extension of the writer’s 3sdimensional 
analysis. However, Rosen's analysis leads him to reaffirm 
the conclusion that one should retain the linear speed- 
distance law, despite the fact that this limits the size of a 
disk which can rotate with a given angular velocity. This 
limitation arises, of course, from the fundamental condi- 
tion that the parts of the disk cannot exceed the speed of 
light, and is in no way a result of the physical properties 
of the material. Since I find myself still unable to accept 
this conclusion, it seems pertinent to consider the question 
again briefly. 

It seems to me that the basic problem here is not howa 
given physical disk would behave if set in rotation. The 
more fundamental question is whether one can find a group 
of transformations representing relatively rotating co 
ordinate systems; in the sense of the theory of continuous 
groups any member of the set should then represent a ref- 
erence system in uniform relative rotation with respect to 
any other member of the set. This problem has an obvious 
solution in classical kinematical theory, but is unsolved 
as yet in relativity theory. The classical group must appear 
naturally as a degenerate case for an infinite value of the 
speed of light. 

Starting with an extended body of any size in a given 
“laboratory” system, on passing to another member of the 
family of reference systems the body should appear to be 
in uniform rotation in the sense in which such motion was 
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defined by the transformation group. It is difficult to imag- 
ine how any condition could arise to limit the size of the 
body in such a treatment if the group possessed infinitesi- 
mal elements, and so could describe rotations with in- 
finitesimal angular velocity. 

Although this version of the rotation problem underlay 
my earlier attack, I have not yet been able to give it an 
acceptable solution, and may take this occasion to pose it 
as one needing an answer. 


1N. Rosen, Phys. Rev. 71, 54 (1947). 
2 EL. Hill, Phys. Rev. 69, 488 (1 946). 


Remarks on C. Frondel’s Letter: “Elastic Defi- 
ciency and Color of Natural Smoky Quartz” 
K. PRrzIBRAM 


11. Physical Institute, The University, Vienna, Austria 
January 27, 1947 


AVING just seen C. Frondel’s interesting letter! I 
should like to make two remarks. 

1. Anomalies after irradiation with radium and x-rays 
have also been observed in the piezoelectric constant of 
quartz, this constant in some cases being increased by 
irradiation.? 

2. The fact that the natural color of smoky quartz is 
more stable than the color produced artificially in irradi- 
ated quartz does not contradict the theory that this natu- 
ral color is of radioactive origin. 

As I have frequently pointed out, orily the most stable 
color and luminescence centers will be produced in nature 
owing to the very small intensity of the radiation acting 
during long periods, as compared with the strong artificial 
irradiation during a short time. I have called this “the 
principle of natural selection of the most stable.’’* 

That colorless or relatively light colored bands in natu- 
ral smoky quartz are more deeply affected by artificial ir- 
radiation shows that they contain more unstable centers 
than stable ones. 


iC. Frondel, Phys. Rev. 69, 543 (1946). 
2 J. Laimbick, Wien. Anz. (May 18, 1928); F. Seidl, Wien. Ber 

142, 467 (1933); F. Seidl and E. Huber, Zeits. f. Physik 97, 671 1985)" 
1K, Przibram, Zeits. f. Physik 68, 413 (1931). 


The Magnetic Quenching of Superconductivity 


M. J. Stenxo AND R. A. OcG, Jr. 
Depariment of Chemistry, Stanford University, California 
February 1, 1947 


HE only previously expressed generalization regard- 

ing temperature dependence of magnetic threshold 
fields for super conductors’? appears to be the often- 
quoted statement that the functional relation is roughly 
parabolic for a given substance. It has been noted that the 
curves for mechanically “‘soft’’ elements (as well as for a 
few alloys) tend to display smaller initial slopes than those 
for “hard” elements (and most alloys) of similar upper 
critical temperature (vide infra). More careful examination 
of the threshold curves for the “soft” superconductors 
Pb, Hg, Sn, In, Tl, CuS, Au,Bi, Zn, and Cd indicates that 
each is reproduced with remarkable accuracy by a law of 
the form Hr=A(T¢!—T!). Here Tc represents the upper 
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critical temperature—i.e., at zero field, whereas T' repre- 
sents the critical temperature in a longitudinal field of Hr 
gauss. A remarkable generality is indicated by the fact 
that the parameter A is found to have nearly the same 
value for all of the substances listed, the deviations being 
of the same order as those for the individual curves. For 
the initial slope at T= T¢, one derives 


(dHr/dT)re= —§AT oA. 


The practical identity of A for the various substances is 
most readily verified by use of this corollary law, since the 
initial slopes are generally measured with greatest accuracy. 
For the listed substances (with a range of T¢ from 7.26 
to 0.54°K) the average value of A is some 50 gauss- 
degree“; with a mean deviation of only a few percent. 
Studies of thin films and wire, as well as of colloidal 
particles'? of “soft” superconductors, have shown the 
threshold fields to increase for sufficiently small linear 
dimensions, by a factor which is the greater, the closer the 
approach to the critical temperature. This suggests that 
mechanical sufficient to cause a micro-lamellar, fibrous, or 
mosaic structure in a massive sample should result in 
serious deviations from the above law—in the sense of an 
increase in the parameter A, differentially greater the 
higher the temperature. The threshold curves for the few 
“hard” superconductors subjected to careful study seem 
to bear some such relation to those of “soft” supercon- 
ductors of similar upper critical temperature, the ratio of 
initial slopes being much greater than the ratio of extrapo- 
lated threshold fields at the absolute zero. It is proposed 
that the above law represents the limiting “‘ideal’’ behavior 
for a truly homogeneous specimen of macroscopic di- 
mensions. The apparent exceptions displayed by the more 
numerous class of “hard” superconductors are attributed 
to irreducible structural inhomogeneities on a micro-scale. 
The incompleteness of the Meissner effect (expulsion of 
magnetic flux) in typical “hard” superconductors is further 
evidence of such structural inhomogeneities. 

A law of such simplicity as the above suggests that the 
electronic nature of the superconducting state may prove 
to be simpler than has usually been proposed. In par- 
ticular, a satisfactory model must presumably allow the 
experimental value of the parameter A to be found by com- 
bination of universal constants only. 


1E. F. Burton, H. Grayson Smith, and J. O. Wilhelm, Phenomena at 
the Temperature of Liquid Helium (Reinhold Publishing Corporation, 
New York, 1940), 
Shoenberg, Superconductivity (Cambridge University Press, 


Decay of Mesons Stopped in Light Materials 


T. SIGURGEIRSSON AND A. YAMAKAWA 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 
February 8, 1947 


E are carrying out an experiment to observe the 

decay of mesons stopped in different light materials. 
The preliminary results seem to indicate’ that there are 
more decay electrons emitted for each stopped meson if 
the meson is brought to rest in material of low atomic 
number than if it is stopped in aluminum or some material 
with a higher atomic number. 
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ABSORBER 


Fic, 1. Experimental arrangement. 


The experimental arrangement is shown in Fig. 1. The 
moderator is The G-M counters for de- 
tecting the mesons and the decay electrons all have 1” 
diameter and 12” sensitive length. The mesons entering the 
absorber are observed as a double coincidence in counter 
trays I and II. The decay electrons coming out from the 
absorber are detected in counter tray III containing four 
counters. 

Absorbers of the following materials have been used: 
Be, C, NaOH, Al, SiC, and S. The decay curves for the 
different absorbers are still rather inaccurate and do not 
show any definite disagreement with a mean life of 2.2 
usec. 

To compare the number of decay electrons from the 
different absorbers, we count the number of impulses 
from counter tray III that are delayed between 1 and 6 
usec. with respect to the coincidences I XII. The most re- 
liable results come from a comparison of beryllium and 
sulphur absorbers of equal stopping power. The same 
number of mesons will be stopped in both, and in both 
cases the same fraction of the decay electrons will be 
stopped in the absorber before reaching the counters. We 
have counted the delay pulses with and without absorber, 
changing absorber every day. An 8-day long run gives 
the result of Table I. 

Similar runs have been made for other absorbers. The 


TABLE I. 
Absorber Be s None 

Mass 3290 g 3440 
Relative stopping power per g/cm* 1.05 (1.00 
Relative number of mesons calcu- 

lated to be stopped 1.00 (1.00) 
Delayed counts 219 213 29 
Time of observation 58 hr. 84.3 hr. 44 hr. 


Counts per hour corrected for 
background 3.1240.18 1.87+0.10 


Relative number of disintegrations 


for  - number of mesons stopped 1.67 +0.14 
Ratio of total number of mesons at 

sea level to number of positively 

charged mesons* 1.8 


* D. J. Hughes, Phys. Rev. 57, 592 (1940). 


TABLE II. 


Absorber Al Cc NaOH sic 
4680g  2730g 3440 & 378, 


Mass 
Counts, corrected for back- 
ground, and reduced 
to equal numbers of 


mesons stopped (1.00) 1.7202 1.4401 1.040, 


results are given in Table II. To get the ratios for the 
number of emitted decay-electrons the figures of Table I] 
have to be corrected for the absorption of the decay. 
electrons. The detailed evaluation of this effect has not yet 
been completed, but the ratio will have to be reduced by 
an amount estimated to be roughly of the order of 20 per- 
cent for C and 10 percent for NaOH and SiC. 

The present results confirm the findings of Convers, 
Pancini, and Piccioni,' that negative mesons absorbed in 
carbon emit decay-electrons whereas no decay-electrons are 
observed if they are absorbed in iron. The absence of decay- 
electrons is easily explained by assuming that the nega- 
tive mesons are captured by the nucleus.? 

Our results indicate that the capture probability (a) in- 
creases gradually with increasing atomic number, but 
(6) disagrees completely in absolute value with the pre- 
dictions’ of the meson theory of nuclear forces. 

This work has been supported by Navy Bureau of 
Ordnance Contract 7920 with Princeton University. We are 
indebted to the Radio Corporation of America, Princeton 
Laboratories, for development of a special time measuring 
and recording device used in these experiments, which are 
continuing and will be further reported. 


1M. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 71, 209 (1947), 
2 S. Tomonaga and G. Araki, Phys. Rev. 58, 90 (1940). 


Mechanism of Capture of Slow Mesons' 


Joun A. WHEELER 


Palmer Physical Laboratory, Princeton University, 
Princeton, New Jersey . 


February 8, 1947 


HE Rome group* observe that negative mesons 

stopped in carbon undergo radioactive decay but 
that those stopped in iron do not. The measurements of 
Sigurgeirsson and Yamakawa' indicate (1) that for other 
light elements the probability of nuclear capture is like- 
wise small, relative to the probability of radioactive decay 
(1/ro=1/2.15X10-* sec.); and (2) that the ratio of the 
two transition probabilities rises with atomic number, 
passing through the value unity for an atomic number, 
Zo, probably in the neighborhood of Zo~10. 

It follows from these observations that the time re- 
quired for nuclear capture in light elements is 10~* sec. or 
more, far longer than the time taken by the meson to reach 
its lowest Bohr orbit about the nucleus via radiation and 
Auger effect.‘ Consequently it is only for a negative meson 
moving in a K-orbit of a nucleus of charge Zo that nuclear 
capture and normal radioactive decay must be considered 
to have the same probability. For mesons in higher orbits 
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or in the free state direct capture by the nucleus must be 
concluded to have negligible probability (¢=2X10~" cm* 
ev/(energy of meson)). 

In the K-orbit of another nucleus of charge Z the prob- 
ability of nuclear capture may be expected to be greater 
(a) by the factor (Z/Zo) for the larger number of protons 
presumably capable of being transformed into neutrons 
and (b) by the volume concentration factor (Z/Zo)* based 
on the expression (A*/uZe*) for the effective radius of the 
Bohr orbit. Combination of factors gives for the expected 
decay constant for negative mesons stopped in matter, 
1/r=(1/ro)[1+(Z/Zo)*], and for the fraction which send 
out decay electrons the expression, 


An appreciable dependence of the mechanism of dissipa- 
tion of energy upon atomic number will be expected to 
alter the value of the exponent, s=4, deduced above. 

A significant fraction of the moderated negative mesons 
will arrive at the 2s orbit and might be expected to remain 
in a metastable state with a decay probability lower than 
that calculated above. However, the finite extension of the 
nuclear charge inverts the normal order of 2s and 2p terms 
(by ~3000 ev in Al) and gives the 2s state a short life 
with respect to radiation (~2X10~* sec. for Al). Also 
direct transitions from 2s to 1s may occur via transfer 
of energy to a 1s electron (4X10-" sec. in Al). Conse- 
quently only nuclear capture from the K-orbit need be con- 
sidered. 

The Auger electrons given out in 2s—>ls and 2p-+1s 
transitions (0.7 Mev in argon, 1.4 Mev in iron) offer an 
independent means to determine the mass of the meson. 
For nuclei of charge Z >23 pair production is also possible 
in the 2s—>1s transition but has negligible probability. 

For a nucleus provided with 100 Mev of excitation by 
meson absorption, evaporation of 6-15 nucleons (nuclear 
“star” ) will be expected to have a probability ~10* times 
greater than the probability of electromagnetic radiation. 
Stars at mountain altitudes are too frequent for all to be 
accounted for by meson absorption. Whether conversely 
all negative mesons stopped at sea level (~400 observed 
per m* of air per day) can be assumed to produce stars 
(~300/m* day as deduced from photographic plates) 
appears an open question in view of existing experimental 
uncertainties. Neither experimentally nor theoretically can 
one exclude the possibility of some other process for dissi- 
pation of the mass energy of the meson. 

Straightforward application of the law of microscopic 
reversibility gives 10-*°(Z/Z»)*-* cm? as upper limit for the 
(y, u) reaction produced by y-rays in the 100-Mev range, 
and a correspondingly low probability for production of 
mesons by heavy particle bombardment. 

A more detailed article is being prepared. 

came title, Bull. Am. Phys. Soc. 22, No. 1, 

2M. Conversi, E. Pancini, and O. Piccioni, Phys. Rev. 71, 209 


(1947). The author is indebted to Dr. Piccioni for preliminary com- 
of these results. rae 


and A. Yamakawa, preceding letter. 
‘E. Fermi, E. Teller, and V. Weisskopf, this e, give an explicit 
formula for the rate of loss of energy by Auger transitions between orbits 


large quantum number. The author is ted to these authors for 
the opportunity to see their letter before publication. 
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Gamma- and Beta-Ray Energies of Some Radio- 
active Isotopes as Measured by a Thin Mag- 
netic Lens Beta-Ray Spectrometer* 

RALL 
University of Chicago, Chicago, Illinois 
AND RoGerR G. WILKINSON 


Indiana University, Bloomington, Indiana 
January 28, 1947 


THIN magnetic lens beta-ray spectrometer, similar 
to that developed by Deutsch! and his co-workers, 
was constructed by the authors for the purpose of measur- 
ing radioactive isotopes produced at the Clinton pile. The 
instrument was designed primarily for the study of beta- 
radiations between the energies of 0.1 to 3.0 Mev. The 
momentum of the focused electrons is determined by the 
relation P=kI, where J is the current through the lens 
coil and & is a constant depending on the geometry of the 
coil. Since the electron paths are quite complicated, no 
attempt is made to use the spectrometer as an absolute 
instrument; i.e., & is determined experimentally. To de- 
termine k, the F conversion line of thorium B (Hp= 1385) 
and the photoelectron line due to annihilation radiation 
from positron emitters (0.511 Mev) are used. Agreement 
between the two calibrations is obtained to about 0.1 
percent. Half-intensity widths between 2 and 4 percent are 
obtained for conversion lines, and between 4 and 6 percent 
for photoelectron lines. Beta- and gamma-sources as weak 
as 0.005 and 0.1 mc, respectively, have been successfully 
used. 
Table I lists the results of most of the investigations. The 
values given for the gamma-ray energies are accurate to 


Taste I. 
Half- Max. beta-ra Gamma-ray Relative 
Element life energy (Mev) ,energy (Mev) intensity 
*Srs? 55d 1.5+.1 
37d {33 
Ag 225d 0.59 +.05 66 4a 
90 47 
1.40 9 
Sb t.57 
60d 1.72 
1m? 30d? 61 
1.2h 18 +.1 
12.5d 1.05 +-.05 t.54 
*Laue 40h 1.45+.1 2 
49 5 
87 10 
1.65 77 
2.3 6 
120d 0.53 +.03 
.22 4 
1.13 37 
1.22 57 
Pb? 3.3h 0.68 +.03 
* Fission product. 


t Internal conversion of the gamma-ray noted. 


about 2 percent. The estimates of the relative intensities 
of the gamma-rays are only approximate and the values 
for the smaller components should be taken as indications 
only. 

* This document is based on wert poteemned for the Manhattan 
Project, and the information covered will appear in Division 
IV of the Manhattan Project Technical Series as a part of the contribu- 


tion of the Clinton tories. 
1 Deutsch, Elliot, and Evans, Rev. Sci. Inst. 15, 178 (1944). 


he 
II 
y- | 
et 
by 
i, 
in 
re 
y- 
a- | 
n- 
ut 
e 
of 
re 
on 
ng 
re 

7). 
ns 
ut 
of 
er 
ay 
he 
or 
ch 
nd 
on 
ar 
ed 
its 


322 LETTERS TO. 


. Ferromagnetism at Very High Frequencies 


I-Magnetic Iron at 200 mc 


M. H. Jounson, G. T. Rapo, AND M. MALOOF 
Naval Research Laboratory, Washington, D. C. 
February 3, 1947 


T is possible to determine the complex permeability, 
“1—tpz, of a magnetic metal by measuring the phase 
velocity and attenuation in a coaxial line whose center 
conductor is the metal under investigation. The fractional 
deviation of the phase velocity from the velocity of light 
in vacuum, §, which is of the order of 0.5 percent, as well 
as the attenuation, a, are essentially fixed by the perme- 
ability and conductivity of the metal. We have succeeded 
in simultaneously measuring 8 and a with sufficient pre- 
cision to determine and pe.! 

Our experiments have been made at 200 mc in a half- 
wave coaxial resonator part of whose center conductor was 
the substance under examination.? The Q of the cavity 
and its resonant frequency* were measured as a function 
of a d.c. polarizing field, Ho, which was everywhere parallel 
to the high frequency magnetic field. In Fig. 1 we have 
plotted our experimental results for magnetic iron. The 
top curve is the ratio of Q at a polarizing field Hy to Q 
= Qbias at a polarizing field of 6.25 oersteds which was con- 
tinuously maintained throughout the measurements. The 
saturation limit, Q./Qpies, was found independently by 
measuring the cavity Q with a brass center conductor and 
making the proper correction for the difference in con- 
ductivity between brass and magnetic iron. The lower 
curve is the difference, 5, in resonant frequency for the 
cavity with a polarizing field Hy and with a polarizing field 
of 6.25 oersteds. Here there is no independent way of fixing 
the saturation limit, 5.., which must be found by extrapola- 


VV 


Fic. 1. Q ratio and change in resonant frequency as a 
function id He. 


THE EDITOR 


tion, Hence, it is much less accurately determined (~42 
percent) than the limit of the upper curve. 
From the ratios Q/Qpiss and Q../Qpiss and the 

of the cavity, a can readily be calculated. Similarly, from 
é and 6. 8 can be calculated. With the known conductivity 
of iron, 4: and yz can then be found. In Fig. 2 we have 
plotted u: and us determined in this way from the Curves 
of Fig. 1. The error in yw; is estimated as less than +5 per- 
cent. However, the values of u2 (which depend on 8) 


are considerably more uncertain. The error in js is esti. 


mated as less than +10 percent for small values of Hy and 
as less than +50 percent at Hy= 500 oersteds. 

In Fig. 2 we have also plotted values of the incremental 
d.c. permeability* for a sample of magnetic iron cut from 
the same rod as the specimen used above. It will be noted 
that the experimental values of u, are in excellent agreement 
with the d.c. incremental permeabilities for Ho>15 oersteds, 
Below 15 oersteds the d.c. permeability is larger than a, 
the ratio reaching a value in excess of 3 at 2 oersteds, 
Values of yw: for Ho<6.25 oersteds, which are not shown in 
Fig. 2, depend very much on the cycle by which the state 
of magnetization is reached as well as on Ho (hysteresis 
effects). However, our measurements show that 4; is al- 
ways less than 60, so that wu: at Ho=6.25 oersteds is within 
20 percent of its maximum value. The measurements also 
show that u2/m: is always less than }. 

These results can be interpreted by a modification in the 
standard description of magnetization: reversible and 
irreversible displacements of the boundary walls between 
domains below the “knee’’ of the magnetization curve and 
rotation of the domain spins toward the direction of the 
field above the knee. Our results show that p2/(u:—1) 
~0.15 and that mw: is very nearly equal to the d.c. incre 
mental permeability above the knee (H»>15 oersteds), 
We conclude that at 200 mc magnetization by rotation of the 
domain spins takes place with a relaxation time, r, approxi- 
mately equal to 10~ sec. (wr~0.15). Our measurements 
indicate that r does not change by more than 50 percent 
in the range 15 <H»)<500 oersteds. 

For small values of Ho the incremental permeability is 
much larger than Since pu2/ui<0.5, the discrepancy 
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Fic. 2. Experimental permeabilities as a function of the 
polarizing field He. 
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cannot be explained by a longer relaxation time in these 
fields. (With »2=0.5 w: the maximum discrepancy could 
be 15 percent instead of the observed factor 3.) The small 
increase in ys below 6.25 oersteds could be entirely due to 
magnetization by rotation of the spin directions. We there- 
fore conclude that at 200 mc magnetization by the displace- 
ment of domain boundary walls is greatly reduced and is 
almost wholly out of phase with the magnetizing field, and 
that magnetization by rotation is the dominant effect. 

This investigation is being extended to other materials 
and to other wave-lengths. 


1 Abstracts L4 and L5, Bull. Am. Phys. Soc., Jan. 30, 1947. 

2 We are indebted to Mr. R. A.C hegwidden of the Bell Telephone 
Laboratories for samples of several magnetic materials. 

#In practice the frequency is held constant and the resonant wave- 
length of the cavity altered linearly by displacing a dielectric bead. 

e are much indebted to Mr. E. A. Gaugler of the Naval Ordnance 

Laboratory for the preliminary measurements of d.c. incremental per- 
meability shown in Fig. 2. 


A Defense of the Cauchy Relations 


CLARENCE ZENER 


Institute for the Study of Metals, The University of Chicago, 
Chicago, Illinois 


February 5, 1947 


oe N' has recently presented a theoretical analysis 
of the Cauchy relations between the coefficients of 
elasticity. From this analysis he concludes that the assump- 


. tion of central forces does not necessarily lead to relations 


of the Cauchy type, and implies that the failure of the 
Cauchy relations in particular lattices has no theoretical 
significance. The purpose of the present letter is to point 
out that the Cauchy relations do follow when in addition 
to the assumption of central forces we assume that each 
atom is at a center of symmetry of the lattice. Since many 
metals and salts satisfy this symmetry condition, namely 
simple cubic, face centered cubic, and body centered cubic 
crystal structures, a failure of the Cauchy relation between 
their coefficients of elasticity is of theoretical significance, 
namely this failure implies that all the forces are not of the 
central type acting along lines joining lattice points. 

The proof that the Cauchy relation follows from the 
above two assumptions has been given by Love.? An al- 
ternative proof is in fact furnished by Epstein’s own analy- 
sis. Epstein shows that the correct relations are identical 
with the Cauchy relations save for the presence of certain 
expression of the type 


= 
where the symbols in the right member are defined in 
reference 1. The above expression may be written in the 
more familiar form 


= 42 pol (Gur), 
where yy and my, are the direction cosines of the vector 
passing from the lattice point « to the lattice point », re- 
ferred to the j and k axes, respectively. Love? has given a 
very simple interpretation to this second expression for 
vjz in the case of central forces when each atom is at a 


_ center of symmetry; namely, he has shown that aside 


from a multiplicative constant, v,z is simply jk, the stress 
acting across a plane normal to the j axis in the direction 
of the & axis. Under these two conditions the Cauchy 
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relations are therefore valid, provided the specimen is 
under no initial stress. 

Mathematical 


5A. heory of Elasticity (Cambridge Uni- 
versity iI 1906), second edition, p. 535. 


The Upper Energy Limit of the K” 
Beta-Ray Spectrum 
W. J. HENDERSON 


National Research Council, Ottawa, Canada 
February 8, 1947 


ONOPINSKI' has discussed the disintegration on 

K* in terms of the selection rules, which are assumed 

to hold during the beta-decay. The nuclear transition in- 

volved is a highly forbidden one, and therefore interesting 

from the theoretical point of view. However, no decision 

could be made between the Gamow-Teller and Fermi 

selection rules because of the conflicting experimental 

values for the upper limit of the spectrum. For this reason 

it is felt that a fuller account of my measurements? of the 
upper limit may be of interest. 

The measurements were made by the absorption method 
wherein the absorption of the K** beta-rays in copper was 
compared with the absorption of radium EZ and uranium 
X. beta-rays under the same experimental conditions. 
Since the activity of potassium is very small, a large solid 
angle and a sensitive detector are required. A cylindrical 
Geiger-Miiller counter of diameter 1.7 cm and length 5 
cm of Dow-metal wall, of 0.0354 g/cm* was surrounded 
by a cylinder of KCI, of inner diameter 3.2 cm and 2.3 cm 
thick. In this arrangement the initial count, with only the 
counter walls absorbing, was 500 per minute. The absorp- 
tion curve obtained by dropping copper tubes over the 
counter is shown in Fig. 1. After the beta-rays have been 
stopped at a range of 0.375 g/cm? of copper, there remains 
a constant background due to the gamma-radiation of 
energy 2.0 10* volts emitted* by potassium, with an in- 
tensity of 3 quanta per 100 disintegrating potassium atoms. 

With the ordinary radioactive bodies, where the in- 
tensity is much greater, it is customary to use some 
empirical relation*® to calculate the upper energy limit 
from the range or the absorption coefficient. Such a pro- 
cedure is not permissible in the case of potassium because 
of the low initial activity, and the range 0.375 g/cm* must 
be corrected by comparison with other radioactive bodies 
where the energy distribution is known. The measurements 
were therefore repeated using identical sized sources of 
radium E and U;O, mixed with NaCl as filler to make up 
the same mass per unit area as the KCI source. In each 
case the strength was adjusted so that the initial intensity 
was approximately the same as the KCl—i.e., 500 counts 
per minute. The radium E and uranium X; curves are in- 
cluded in Fig. 1. From these absorption curves the ap- 
parent ranges of the radium E and uranium X; beta-rays 
are 0.310 g/cm? and 0.865 g/cm*, respectively. The correct 
ranges of the beta-rays from these bodies are very well 
known,** and are 0.475 g/cm? and 1.10 g/cm*. It is ap- 
parent that the ranges measured with low intensity sources 
are too small, and the usual relations cannot be used to 
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find the upper limit. From the discrepancy between the 
observed ranges and the true ranges for radium E and 
uranium X:2, we can estimate that the correct range for 
the potassium beta-rays is 0.55 g/cm’. This value, by the 
use of Feather’s’ range energy equation R=0.543XE 
—.160, gives 1.310 volts as the upper limit of the K* 
beta-ray spectrum. 

No claims of exactness can be made for this value, but 
it is clear from the curves that K* emits beta-particles of 
higher energy than does radium E, where the upper limit 
is known*" to be 1.15 10* volts. Another estimate may 
be made from the absorption coefficients since the absorp- 
tion curves are exponential over a large part of their 
length. The absorption coefficients of radium E and ura- 
nium X:, and data from other bodies where the upper 
limits and absorption coefficients are known give a rough 
u/p vs. upper limit curve. From this curve the u/p=12.8 
for potassium yields an upper limit of 1.4X10* volts. The 
last estimate is valid only if the shape of the potassium 
spectrum is sensibly the same as that of the other-bodies. 
That the K* spectrum is not highly complex is indicated 
by the low intensity of the gamma-rays emitted. It is felt 
that the value derived from the corrected range is more 
reliable and that the best value is 1.3 10* volts with an 
estimated error of 0.15 X 10° volts. 

The former lower values of about 700 kilovolts for the 
upper limit of K*, derived largely from absorption data, 
are due in part to the fact that corrections were not made 
for the low activity of the source. Libby and Lee” have 
investigated the K* spectrum with a concentric arrange- 
ment of source and counter, and an axial magnetic field, 
and give a value of 725 kilovolts, but state that no effort 
was made to obtain an accurate value. This type of meas- 
urement with a more elaborate system of counters, has 
been repeated recently by Dzelepow, Kopjave, and 
Vorobjov," who find a value of 1.35 X 10° volts. 

These measurements were done some years ago at 


% TRANSMISSION 


° o2 08 of os 


Fic. 


THE EDITOR 


Queen’s University, Kingston, Canada, and | wish to 
thank Dr. J. A. Gray, F.R.S., for the use of the equipmen 
and radioactive sources used in the experiment. . 


1 E. Konopinski, Rev. Mod. P’ 15, 209 (1943), 
2 W. J. Henderson, Phys. Rev. 55, 238 (1939). 


40. Klemperer , Proc. Roy. Soc. A148, 638 (1935), 
5N. Feather, Proc. Phil. Soc. 34, 599 (1938), 
* Rutherford, Chadwick, and Ellis, Radiations from Radioactive Sub. 


sonces 
Nt Reather, Proc. Camb. Phil. Soc. 34, 115 (1938), 
oR W. Sargent, Proc. Roy. Soc. A139, 659 (1933), 
fig ant G. Ward, Can. J. Research 15, 42 (1937), 
10 A. Flammersfeld, Physik. Zeits. Lg (1937), 
u E, M. Lyman, Phys. Rev. 51, 1 (193 7). 
2 W. Libby and 2 Lee, Phys. Rev. 55, 245 (1939). 
capenh. Dzelepow, M. Kopjava, and E. Vorobjov, Phys. Rev. 69, 538 


+L. : Gray and G. T. P. Tarrant, Proc. Roy. Soc. A143 695 (1934) 


Coincidence Experiments on Ga” 


ALLAN C. G. MitcHELL, Ep. T. JURNEY, AND MARGARET 
Indiana University, Bloomington, Indiana 
January 29, 1947 


A PRELIMINARY investigation has been made, with 
coincidence methods, of the radiations from the 14,9. 
hour radioactive Ga™. A sample of gallium was activated 
by exposure to slow neutrons from the cyclotron. The 
short period activity (20 min.) of Ga’ was allowed to die 
out and the longer period activity investigated. The 
radiation consists of both gamma- and beta-rays. 


The most energetic gamma-ray was found by measuring | 


the energy of Compton recoil electrons ejected from an 
aluminum radiator. The recoil electrons were made to 
pass through two thin-walled counters connected in a co- 
incidence circuit. Thin sheets of aluminum were placed 
between the counters until the number of coincidences 
dropped to a small, constant background. The range of the 
recoil electrons was found to be 1.22 g/cm® (including 
counter wall thickness). From the curve of Curran, Dee, 
and Petrzilka! for the energy of the gamma-ray in terms 
of the range of the Compton electrons ejected from alumi- 
num, the energy of the most energetic gamma-ray is found 
to be 2.4+Mev. An absorption curve of the gamma-rays 
in lead shows, in addition to the hard component, a soft 
component whose absorption coefficient corresponds to 
that of a gamma-ray of about 0.82 Mev energy. 

The beta-ray endpoint was determined by absorption 
in aluminum and was found to be 2.3 Mev. 


° 4 & 3 4 5 & 7 8 
Total Absorber gm/cm*® 
Fic. 1. 
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Experiments were performed on gamma-gamma- and 
beta-gamma-coincidences.* The source was placed between 
two counters in coincidence. Enough aluminum was placed 
on each side of the source to stop all beta-rays. Gamma- 
gamma-coincidences were measured and the number of 
coincidences per recorded gamma-ray was found to be 
0.39 10-*. Aluminum was now removed from one side of 
the source and the number of beta-gamma-coincidences 
measured as a function of the energy of the beta-rays. The 
results are shown in Fig. 1, in which the number of beta- 

ma-coincidences per recorded beta-ray is plotted as a 
function of the thickness of absorber in g/cm’. It will be 
seen at once that the beta-ray spectra of Ga” is complex. 

Since beta-gamma-coincidences were observable prac- 
tically out to the beta-ray endpoint, it follows that the 
beta-transition does not lead to the ground state of the 
product Ge”. In addition, there is a softer group of beta- 
rays having an endpoint at approximately 0.77 Mev. 

From the evidence available—beta-rays of 2.3 and 0.77 
Mev; gamma-rays of 0.8 and 2.4 Mev—no reasonable 
energy level scheme can be proposed unless one assumes 
that the low energy gamma-radiation consists of two 
quanta. Miller and Curtiss* have found three gamma-rays 
0.64, 0.84, and 2.25 Mev. Assuming these three gamma- 
rays and using the results of the coincidence experiments, 
the following level scheme is possible. The 2.4-Mev beta- 
ray leads to an excited state of Ge” which is 1.46 Mev 
above the ground state. The 0.64- and 0.82-Mev gamma- 
rays are emitted in cascade, from this level. The 0.77-Mev 
beta-ray goes to a level 3.3 Mev above the ground level, 
from which the 2.4-Mev gamma-ray leads to the 0.84-Mev 
level and thence to the ground state. 

This work was carried on under contract with the Office 
of Naval Research. 


1Curran, Dee, and Petrzilka, Proc. Roy. Soc. 169, 269 (1938). 

2For a detailed description of this experimental technique, see 
, Mitchell, and McDaniel, Phys. Rev. 56, 422 (1939). 

3 Miller and Curtiss, PPR 9B, 7 -45 -6 (1946). 


Non-Primary Cosmic-Ray Electrons above 
the Earth’s Atmosphere 


G. J. PERLow J. D. SHipman, Jr. 
U. S. Naval Research Laboratory, Washington, D.C. 
February 7, 1947 


N a V-2 rocket fired to a height of 70 miles on January 
10 at White Sands, New Mexico, data were obtained 
on the penetration of cosmic rays through 2 and 4 cm of 
lead. In addition showers below each lead slab could be 
measured. The apparatus is diagrammed. schematically in 
Fig. 1. Its axis was at 45° to that of the rocket and during 
flight it.precessed about the latter. The rocket axis: was a 
few degrees off vertical. The telescope pointed through a 
steel “window” 2.8 mm thick. Protection against the war- 
head showers, found in this and each previous flight, was 
obtained by anti-coincidence counters (marked with crosses 
in the figure). The arrangement of the equipment in the 
warhead was such that a considerable amount of lead 
shielded most of the solid angle from below. The signifi- 
cance of this will be seen presently. 
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While a full account of the experiment awaits more de- 
tailed analysis, it seems desirable to report one interesting 
result at this time—namely, the presence above the earth's 
atmosphere (<1 mm Hg) of relatively large numbers of 
particles absorbed in the lead. About 25 percent of the total 
ionizing component was stopped in the first 2 cm and 10 
percent in the second. Of the latter group about 1 in 8 
produced showers tripping at least 3 counters under the 
first slab. The number of particles stopping in 4 cm or 
less was about the same as the number which penetrated 
without multiplication. The remainder of the radiation 
penetrated 4 cm and produced showers either under 2 or 
4 cm, or under both. The soft particles (range =4 cm) had 
an intensity in free space roughly 15 percent of the total 
intensity observed in the flight at the Pfotzer maximum. 

It seems most reasonable to ascribe this soft radiation 
to electrons. On the basis of the, cascade theory these will 
have energies <5 X 10* ev for the most part. An extrapola- 
tion by Mr. Siry of this laboratory of the Lemaitre-Vallarta 
curves to 40°N geomagnetic latitude gives as the minimum 
energy for an extra-terrestrial electron ~4-5X 10° ev at 
45° zenith angle. Thus the soft electrons appear to come 
as re-entrant particles generated in some atmospheric 
layer below. 

The existence of such a component was suggested to the 
group at this laboratory by Professor J. A. Wheeler in a 
private communication some time ago. It was pointed out 
that the almost isotropic angular distribution of the decay 
electrons from low energy mesons would result in some 
vertical intensity upwards and that these would describe 
helical oribits about the magnetic field lines. It is possible 
that the soft particles observed are to be ascribed to this 
hypothesis. ; 

If the downcoming soft particles observed arise from 
some atmospheric layer below, then an appreciable upward 
intensity must exist and it would appear that the results 
of certain types of balloon experiments reported in the 


Fic. 1. Counter telescope in V-2 warhead. 
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literature might require some reinterpretation. A search TasLe I, Ammonia inversion spectrum. 
will be made in a future rocket for an upward intensity. =—- 
Observed Calculated 
E. Golian, E. H. Krause, and G. J. Perlow, Phys. Rev. 70,223 J  K frequency fi Culated minus 
(1946), and 70, 176 (1546). Tequency observed frequency 
138.83 
M. W. P. STRANDBERG, R. T. WENTINK, JR., AND R. E. HILLGER 28 90 
Research Labor of Electronics, Massachusetts Institute of 2 1 23,098.78 23,098.87 re s 
Tec gy, Cambridge, Massachusetts 3 2 22,834.02 22,834.17 +0.15 
February 8, 1947 4 3 22,683.73 22,686.40 +2.67 
5 4 22,653.00 22,652.81 —0.19 
RECISION measurements of frequencies of some of § 
the lines of the ammonia inversion spectrum in the = 33,232.07 —0.09 ‘ 
microwave region have been made. Previous measure- 10 9 24,205.20 24,206.93 py — 
ments'~* have been limited in accuracy to two to five 3 3223449 
mc/s. By the use of a 50 kc/s standard crystal, we believe 5 3 21,285.32 21,291.75 +6.43 Lo 
6 4 20,994.63 20,993.91 —0.72 
that the present measurements are good to at least +50 7 5 20,804.76 20,804.14 —0.62 
kc/s. Because of their intensity these lines may well be = 
* useful as secondary frequency standards. 20 50 20,851.60 —0.90 
There has also been considerable interest'“* in obtaining 42 49  21'301.55 = 
a function to describe these frequencies in terms of the $ 1 
rotational angular momenta. We have calculated the con- 6 3 19,757.55 19,735.12 2243 
7 4 19,218.52 19,217.73 -0.79 


stants in such a function with parameters, J, the total 
angular momentum quantum number, and K, the angular ——— 
momentum quantum number referring to the symmetry 

axis, in the form suggested by the paper‘ of H. Y. Sheng, ”~”° exp(y/ vo) +6, 

E. F. Barker, and D. M. Dennison. The form of this func- v= frequency in mc/s, 


tion which we have found convenient for calculation is as ”°= 23,785.75 mc/s, 
y = —151.450(J?+ J —K*)+59.892K? mc/s, 


follows: 
6=0.01353(J?+J — K*)?+-0.00461 K*(J?+ J —K?) 
| —0.00986K* me/s, L 
= aihiee The data are given in terms of frequency since this is 
ou*Ke2 . our fundamental unit of measurement. The measured and slo 
q calculated frequencies are given in Table I. pre 
The deviations from the calculated curve are plotted in 20 
. Fig. 1, excluding the data for the K=3 lines, since these ’ 
A ; lines have apparently an anomalous behavior. These de- oni 
viation curves are quite regular except for K=3 lines. 
: Our formula for the frequency of these lines allows us to pal 
x predict quite accurately, making use of the deviation tro 
curves, the anomalous deviations that appear on the K =3 nec 
lines, and the values are tabulated against the total angular the 
a momentum quantum number in Table II. It should be (wi 
5.6 gg 8 9 0 K TABLE II. Anomalous frequency shift in K =3 lines. Pu 
the 
“0. J K vmeasured—¥predicted one 
Z 3 42139 EX) 
“1 noted that there is no apparent shift of the K=6 and K=9 am 
Ts . lines, from a fairly smooth curve. There are no K=0 — 
transitions in this inversion spectrum. a 
8 * This research has been supported in part through the Joint Service at 
Contract No. W sc-32037. 
- 1B. Bleany and R. P. Penrose, Nature 157, 339 (1946). 
? W. E. Good, Phys. Rev. 70, 213 ver 
*C. H. Townes, Phys. Rev. ity 
4H. Y. Sheng, E. F. Barker, D. M. Dennison, Phys. Rev. 60, 786 Lat 


Fic. 1. (1941). 
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